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ABSTRACT 

We present deep near-infrared J^, H, and band imaging of a field around MS 1054-03, a massive cluster 
at z = 0.83. The observations were carried out with the ISAAC instrument at the ESO Very Large Telescope 
(VLT) as part of the Faint InfraRed Extragalactic Survey (FIRES). The total integration time amounts to 25.9 h 
in Js, 24.4 h in H, and 26.5 h in band, divided nearly equally between four pointings covering 5'.5 x 5'.3. 
The 3(7 total limiting AB magnitudes for point sources from the shallowest to the deepest pointing are /'°ab = 
26.0-26.2, H^^ = 25.5-25.8, and /T^g = 25.3-25.7 mag. The effective spatial resolution of the coadded 
images has FWHM = 0."48, 0"46, and 0."52 in 7s, H, and K^, respectively. We complemented the ISAAC data 
with deep optical imaging using existing Hubble Space Telescope WFPC2 mosaics through the F606W and 
F814W filters and additional U, B, and V band observations we obtained with the VLT FORSl instrument. We 
constructed a ^Tj-band limited multicolour source catalogue down to Kl°^g « 25 mag 5(t for point sources). 
The catalogue contains 1858 objects, of which 1663 have eight-band photometry. We describe the observations, 
data reduction, source detection, and photometric measurements method. We also present the number counts, 
colour distributions, and photometric redshifts Zph of the catalogue sources. We find that our A's-band counts 
at the faint end 22 < K^^ab S 25, with slope dlogN/dm = 0.20, lie at the flatter end of published counts in 
other deep fields and are consistent with those we derived previously in the Hubble Deep Field South (HDF- 
S), the other FIRES field. Spectroscopic redshifts Zsp are available for w 330 sources in the MS 1054-03 
field; comparison between the Zph and Zsp show very good agreement, with (|zsp-Zph|/(l +Zsp)) = 0.078. The 
MS 1054-03 field observations complement our HDF-S data set with nearly five times larger area at « 0.7 
brighter magnitudes, providing more robust statistics for the slightly brighter source populations. 
Subject headings: cosmology: observations — galaxies: evolution — infrared: galaxies 



1. INTRODUCTION 

Deep imaging surveys for high redshift studies have, in 
recent years, been carried out over an increasing range of 
wavelengths from the X-ray to the radio regimes. An impor- 
tant motivation is the recognition that surveys at any given 
wavelength range may provide an incomplete and biased cen- 
sus of high redshift galaxy populations. Moreover, a multi- 
wavelength approach is highly desirable to assess the nature 
of and establish links between high redshift samples selected 
at various wavelengths, and to constrain the properties of 
their stars and gas contents as indicators of their evolution- 
ary state. Deep surveys at near-infrared wavelengths (NIR; 
A = 1-2.5 /im) have emerged as an important way to investi- 
gate the high redshift universe (see, e.g., McCarthv 2004 for 
a review, and references therein). NIR observations allow one 
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to access the rest-frame optical light of galaxies at z 1-4. 
NIR surveys provide a vital complement, for instance, to op- 
tical surveys, which probe the rest-frame UV light at these 
redshifts fe.g. lShaplev et al. l200UlPapovich et al.l2001l) : flie 
rest-frame optical is less affected by the contribution from 
young massive stars and by dust obscuration, and is expected 
to be dominated by longer-lived lower-mass stars. Selection 
in the NIR, compared to the optical, thus provides a more reli- 
able means of tracing the bulk of the stellar mass at high red- 
shift, of probing earlier star formation epochs, and of identify- 
ing galaxies with spectral energy distributions (SEDs) similar 
to those of galaxies locally. A significant advantage of NIR 
surveys is that they enable, in the same rest-frame optical, di- 
rect and consistent comparisons between z ~ 1 - 4 samples 
and the well-studied populations at z < 1 . 

With the sensitive instrumentation and large format NIR 
detectors mounted on ground-based 8 -10 m telescopes or 
space-based facilities available, it is possible to carry out 
NIR surveys of comparable depth and quality as in the op- 
tical. This greatly helps in the identification of counterparts 
in optical observations and in the construction of accurate 
SEDs. Yet, NIR surveys with sufficient depth and angular 
resolution remain challenging. For instance, a typical nor- 
mal nearby L* galaxy would have an appare nt A'-band mag- 
nitude ^Tveea ~ 23 mag if placed at z ~ 3 ( Kuchinski et al. I 
1200 Ih . Among current NIR surveys, those reaching such faint 
limits are still scarce. They include the Hubble Space Tele- 
scope (HST) NICMOS fiel ds of t he Hubble Deep Field North 
(HDF-N;|Thompson et al. i ll999tj Dickinson et al. 2003) and 
Hubble Ultra Deep Field (UDF: HTiompson et al. ..2005i) . the 
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Hubble Deep Field South (HDF-S) parallel NICMOS field 
((Williams et al. 200§ and the two UDF parallel NICMOS 
fields ( Bou wens et al. Ground-based surveys include 

the Subaru Deep Field (Maihara et al. 2001), a nd several 
blank fields imaged at the Keck Telescopes (Diorgovsk i et al. I 
fr995: Moustakas et al. 1997: Be rshady et al. 1998). 

In this context, we initiated the Fain t InfraRed Extragalac- 
tic Survey (FIRES: iF^nx et al. Il2000h carried out at the ESO 
Very Large Telescope (VLT). The survey consists of very deep 
imaging in the 7s, H , and bands with the ISAAC instru- 
ment ( Moorwood_et al. 1998) of two fields with existing deep 
optical HST WFPC2 data: the HDF-S and MS 1054-03, a 
field around a massive cluster at z = 0.83. In total, 103 h of 
ISAAC observations were spent in a single 2. '5 x 2. '5 pointing 
for HDF-S and 77hina5f5x5f3 mosaic of four pointings for 
MS 1054-03. The K,-hmd 3a limiting total AB magnitudes 
for point sources are « 26.2 mag for HDF-S and w 25.5 mag 
for the MS 1054-03 field, and the effective resolution of the 
coaddedmaps inJ,, H, andK, is 0."45-0"52 FWHM for both 
fields. Together with the HST WFPC2 data and additional 
VLT FORSl optical imaging for the MS 1054-03 field, this 
provides a unique data set with seven or eight band photom- 
etry covering 0.3-2.2 /im over a total area of w 35 arcmin^. 
With the HDF-S data, FIRES reaches fainter limits in all three 
NIR bands than any other existing ground-based survey and 
provides the deepest TiT-band map to date. The MS 1054-03 
field is an essential complement as it considerably increases 
the survey area to an « 0.7 mag brighter limit. 

The data and source ca talogues of the HDF-S are described 
by iLabb e et al. ' ("2003a"). Several key results of both fields 
have been published, including SED modeling, spectroscopy, 
and clustering analysi s of the new popula tion of "Distant 
Red Galaxies" at z > 2 (fFranx et al. 2003; van Dokkum et al. ' 
1^003, 2004; Daddi et al. 2003; Forster Schreiber et al. 2004; 
|Labbe et al. 2005), determination of the stellar rest-frame op- 
tical luminosity density, colours, and mass density at z = 0-3 
dRudnick et al. n 2001 . 2003), determinat ion of the rest-frame 
optical sizes of galaxies out to z ~ 3 (iTruiilloet al.ll20()l 
|2005), and the disc overy of six large disk-Uke galaxies at 
z= 1.4-3.0 (Labbe et al. "2003b'). In this paper, we present 
the data and source catalogues of the MS 1054-03 field. As 
for HDF-S, the MS 1054-03 raw and reduced data sets are 
available from the FIRES Web site'". Sections |2l and |3] de- 
scribe the observations and data reduction. Sect 0] discusses 
the final reduced images. Section|5]describes the /iQ-band lim- 
ited source detection and photometric measurements methods 
and §|6lthe photometric redshift determinations. Section 
lists the catalogue parameters. Section |8] presents the com- 
pleteness analysis, number counts, magnitude and colour dis- 
tributions of the A's-band selected sources. Section |9l sum- 
marizes the paper. All rnagnitu des are expressed in the AB 
photometric system dOkel 19711) unless explicitely stated oth- 
erwise. 

2. OBSERVATIONS 

2. 1 . The MS 1054 - 03 Field 

Our strategy for the FIRES survey was to complement our 
very deep NIR imaging of HDF-S with a wider field, al- 
lowing for better statistics on source populations down to 
flux levels twice brighter than the faint limits of HDF-S. The 
MS 1054-03 field was ideally suited for this purpose because 

' " http ://w w w. strw. leidenuniv. nl/~fires 



of available HST WFPC2 imaging offering a unique combina- 
tion of field size and depth. The WFPC2 data consist of a mo- 
saic of six pointings covering « 5' x 5' observed through the 
F606W and F814W filters (effective wavelengths Aes = 6030 
and 8040 A, and effective widths A A = 1500 and 1540 A, re- 
spectively). " We will denote these bandpasses as Vfeoe and 
/814. The total integration time was 6500 s per band and point- 
ing. To extend the optical wavelength coverage, we obtained 
UBV imaging with the VLT FORSl instrument. 

The MS 1054-03 field benefits from an extensive amount 
of existing data over the electromagnetic spectrum, includ- 
ing a deep mosaic at 850 taken with the Submillimeter 
Common-User Bolometric Array (SCUBA; K. K. Knudsen 
et al. , in preparation), de ep radio 5 Gjizj maging from the 
Very Large Array (VLA; Best et al. 2002), and X-ray ob- 
servations with the ACIS camera onboard the Chandra Ob- 
servatory (|j^^aetar|L2001; Rubin et al. 2004). Optical 
imaging with ACS onboard HST has also now become pub- 
licly available Spectroscopic redshifts have been measured 
for over 400 galaxies in the field from data collected at the 
Keck and Gemini telescop es and at the VLT ( T ran et al. 199^ 
l2?)?)3tlvan Dokkum et al .112 000. 2003, 2004; S. Wuyts et al. , 
M. Kriek et al. , in preparation). These spectroscopic determi- 
nations enabled us to calibrate our photometric redshifts (§|6|l. 

The z = 0.83 MS 1054-03 cluster itself has but a small 
influence over the area surveyed. The majority of galaxies 
detected in the mosaics lie in front of or behind the cluster. 
The cluster mass distribution has been modeled accurately by 
^oekstra. Franx. & Kuiiken (2000) from weak lensing analy- 
sis, so that the gravitational lensing can be well accounted for 
The average background magnification effects over the field 
of view covered by the FIRES observations range from 5% 
to 25% between z = 1 and 4, and are largest in the cluster's 
central regions. We emphasize that colours and SEDs are not 
affected since gravitational lensing is achromatic. Moreover, 
the surface brightness is preserved and, to first order, the sur- 
face density to a given apparent magnitude is insensitive to 
lensing (the exact effect depends on the slope of the luminos- 
ity function of the sources under consideration). 

TableHsummarizes the properties of the MS 1054-03 data 
set and Figure ^ illustrates the fields of view of the differ- 
ent instruments. The observation s and reduction of th e HST 
WFPC2 mosaics are presented bv lvan Dokkum et al. I SlOOd) . 
In this and the following section, we focus on the data taken 
with ISAAC and FORSl, both mounted at the VLT Antu tele- 
scope. All VLT observations were acquired in service mode 
as part of the ESO Large Program 164.O-0612 (PI. M. Franx). 

2.2. ISAAC Near-infrared Observations 

The MS 1054 - 03 field was observed with ISAAC 
dMoorwood et al. l|l998h in imaging mode through the Js, H, 
and filters (Aeff = 1.25, 1.65, and 2.16 fim, and AA = 0.29, 
0.30, and 0.27 /im, respectively). The observations were car- 
ried out on 36 nights, in 1999 December, 2000 May, and 
2001 February through June. The short- wavelength camera 
of ISAAC is equipped with a Rockwell Hawaii 1024 x 1024 
HgCdTe detector array and provides a scale of 0"148 pixel"' 

' ' All bandpass wavelengths and widths quoted in this paper refer to the ef- 
fective transmission profile accounting for the filters and all other instrumen- 
tal optical components, the telescope mirrors, the detectors, and the Earth's 
atmosphere. 

The ACS data have been obtained too recently to be included in the 
work presented here, which has extended over several years and focusses on 
the initial phase of FIRES. 
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for a total field of view of 2. '5 x 2.'5. The WFPC2 mosaic was 
covered with a square of 2 x 2 pointings (referred to as Fl to 
F4; see Figure [Q. Nominal pointings ensured that adjacent 
fields overlap by « 10". 

The data were obtained in "observing blocks" (OBs) con- 
sisting of sequences of 30 exposures totaling 1 h integration 
time in 7s and H, and 40 exposures totaling 40 m in A's. Suc- 
cessive exposures within each OB were dithered randomly in 
a 20" box to allow the construction of sky frames with min- 
imal object contamination. Typical exposure times in J^, H, 
and were of 120, 120, and 60 s, split into 4, 6, and 6 sub- 
integrations, respectively. The entire data set includes 3207 
frames grouped in 93 OBs, with total integration times of 
6-7 h per field and per band. We rejected about 3% of the 
frames, representing 3% of the total integration time, because 
they did not meet our atmospheric seeing or transmission re- 
quirements, or because they were affected by strong instru- 
mental artefacts. Throughout the paper, we consider only the 
3104 frames from 92 OBs with sufficient quality. 

Most of the data were obtained under excellent seeing, with 
median full-width at half maximum FWHM = 0"47 in 7s and 
0"45 in both H and Ks- Ninety percent of the data have a see- 
ing better than 0"65. Figure Ogives the seeing distribution 
of the individual frames. The FWHM was determined by fit- 
ting Moffat profiles to bright, isolated, and unsaturated stars 
in each image (mostly three to five depending on the field, on 
the signal-to-noise ratio of the stars, and on whether they lie 
within the field of view of individual dithered frames), and 
averaging the FWHMs per image. The conditions were pho- 
tometric for about half of the observing nights, or 53% of the 
OBs. 

In the course of the 19-month interval during which the NIR 
observations were taken, the Antu primary mirror was twice 
re-aluminized, in 2000 February and 2001 March. We will 
refer to the three periods separated by these events as periods 
PI, P2, and P3. The effects of re-aluminization are noticeable 
in our data, with an increase in sensitivity between successive 
periods as deduced from bright stars counts in the MS 1054- 
03 data and reflected in the photometric zero point s. B ased 
on the standard stars data used for flux calibration (§ 13. U . the 
increase in 7s, H, and Ks amounts to 37%, 47%, and 50% 
between PI and P2, and 29%, 26%, and 19% between P2 and 
P3, for a net gain of w 80%. Seven OBs were taken in PI, 
41 in P2, and 45 in P3. The ISAAC pixel scale also varied 
slightly over time, by a maximum of 0.27%. 

Figure|3]shows the variations with time and airmass of the 
relative count rates of stars and of the background levels in 
the individual frames. The stars are the same as used for the 
seeing measurements. Their counts were integrated in a 6"- 
diameter circular aperture in the sky-subtracted frames and 
normalized to the median values in period P3. The back- 
ground was computed as the median count rates over the raw 
frames. The figure also plots the seeing FWHM versus air- 
mass. 

The stars count rates increase sharply after each realu- 
minization and, over period P2, are consistent with the long- 
term progressive loss of sensitivity between two realuminiza- 
tions. In contrast, the background levels have overall re- 
mained fairly constant, with no systematic trend with period 
and the variations reflecting mostly the nightly observing con- 
ditions. This constancy in average background levels was also 
noticed in the HDF-S data set and presumably indicates that 
light is scattered in as much as out of the beam. It implies 
a substantial increase in telescope efficiency and in source 



signal-to-noise ratio (S/N) achievable in background-limited 
data after realuminization. 

The NIR background can vary on timescales of minutes; 
it is dominated by telluric OH airglow lines in the 7s and H 
bands, and by thermal emission from the telescope, instru- 
ment, and atmosphere in the Ks band. The nightly variations 
in our data are generally largest and most rapid in the H band 
and show little correlation with the hour of observation. The 
A'-band background has comparable levels but smaller short- 
term variations and is systematically higher at the beginning 
and end of the night. The 7s background is the lowest and 
most stable, and peaks at the beginning of the night. 

Variations in stars count rates within each period reflect 
changes in sky transparency and there is no apparent airmass 
dependence. This indicates that the atmospheric extinction is 
negligible for airmasses of 1.05-2.2 as covered by the ob- 
servations. The seeing is expected to correlate with airmass 
(roughly as FWHM cx [airmass]"'^) but no obvious trend is 
seen in our data, partly due to the small range of airmass and 
because the seeing is largely determined by the actual ob- 
serving conditions. The data suggest an overall increase in 
background levels with airmass, most clearly in H and K^, 
which could be explained by the increasing path length of the 
thermal and OH-line emitting layers. However, for succes- 
sive observations such as within OBs, variations in airmass 
are directly coupled with time so that part of possible trends 
with airmass may be attributed to nightly patterns in the back- 
ground. 



2.3. FORSl Optical Observations 

The MS 1054 - 03 fi eld was observed with FORSl 
(|Appenzeller et al. ' '1998') in direct imaging mode on six 
nights between 2000 January and March using the Bessel 
U, B, and V filters (Aes = 3680, 4330, and 5550 A, and 
AA = 310, 870, and 1070 A, respectively). FORSl employs 
a 2048 X 2048 pixels Tektronix CCD detector. The standard 
resolution collimator was selected for the observations, giv- 
ing a 0"20 pixel"' scale and a 6f8 x 6.'8 field of view cover- 
ing entirely the ISAAC and WFPC2 mosaics. The data were 
collected in OBs of six successive single exposures of 550, 
400, and 300 s in U, B, and V. The total integration times are 
4.6, 2.0, and 0.5 h, respectively, distributed over 5, 3, and 1 
OBs. This excludes the few aborted or complete OBs with in- 
sufficient quality, lying well outside our observing conditions 
specifications. Within each OB, the exposures were dithered 
in a box of w 40" x 30". 

The average seeing FWHM of the data is 0"64 in U, 0"55 
in B, and 0."67 in V. The seeing was determined as for the 
ISAAC data by fitting Moffat profiles to five bright, isolated, 
unsaturated stars in the individual frames. The B- and V-band 
data were all taken prior to the 2000 February realuminization 
while the U -band observations were carried out both before 
and after Our MS 1054-03 data indicate an increase in U- 
band sensitivity of w 20% based on stars count rates in a 6"- 
diameter aperture and small rotation angle differences in the 
camera orientation (0.15° -0.57° between various OBs). The 
data were obtained during photometric nights except for the 
three t/-band OBs taken in 2000 March. 

3. DATA REDUCTION 
3.1. ISAAC Data Reduction 
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We reduced the ISAAC data using the DIMSUM''' pack- 
age within IRAF''*, the ECLIPSE'^ software, and additional 
routines we have developed specifically for our FIRES ob- 
servations. We followed the gener al procedure for the HDF-S 
described by Labbe et al. (2003a) with adjustments necessary 
because of the different characteristics of the MS 1054-03 
data set. In particular, the observations were taken over a 
longer period where the instrument behaviour changed notice- 
ably. Also, the shallower data and more crowded field around 
MS 1054-03 required optimization of various steps to mini- 
mize residual noise and systematic effects. 

3.1.1. Dark, Bias, and Flat Field 

We generated dark+bias images for use in making flat fields 
by averaging together individual dark frames taken with each 
detector integration time (DIT) employed for the flat fields. 
The typical level of the dark frames for a given DIT varied 
significantly over time with five distinct epochs identified in 
our data. We thus created dark+bias images appropriate for 
each epoch in addition to those made from all available dark 
frames. 

We constructed flat fields from sequences of successive ex- 
posures of the sky taken during twilight, on nights where the 
MS 1054- 03 field was observed or as close as possible in 
time, using the ECLIPSE task "flat." This task also produced 
the initial bad (dead) pixel masks. We made three flat fields 
per night and band as follows: (1) by forcing ECLIPSE/flat to 
compute the dark+bias by itself, (2) by providing the average 
dark+bias image from all available data for the corresponding 
DIT, and (3) by providing the dark+bias image for the corre- 
sponding DIT and observing epoch. 

The flat fields often exhibited abrupt jumps at mid-array 
(row 512) with amplitude up to 20% and significant differ- 
ences in large-scale gradients between the three methods, typ- 
ically at the 5% level. These features carry the imprint of the 
detector bias, which dominates the signal in the dark frames. 
The bias level is very uniform along rows but varies impor- 
tantly along columns, with a steep initial decline flattening 
progressively from bottom to top of each of the lower and up- 
per halves of the array frame. The bias depends not only on 
the DIT but also on the detector illumination during and prior 
to an exposure, so that for a given DIT the dark+bias image 
may not represent well the dark+bias present in the twilight 
frames. For each night and band, we chose the flat field for 
which the bias signatures were minimized. This was gener- 
ally the one obtained with method 1 , which better accounted 
for the intrinsic dark+bias of the twilight frames, but in a few 
cases method 3 was the most successful. The nightly flat fields 
also varied noticeably between epochs where the dark+ bias 
changed most but were otherwise very similar. We created 
"master flat fields" by averaging the adopted nightly flat fields 
for each epoch. All flat fields were normalized to unit mean 
over the image. 

3.1.2. Sky Subtraction, Cosmic Rays Rejection, and OB 
Combination 

" DIMSUM is the Deep Infrared Mosaicking Software package developed 
by Peter Eisenliardt, Mark Dickinson, Adam Stanford, and John Wai'd, and is 
available via ftp to ftp://iraf.noao.edU/iraf/contrib/dimsumV2/dimsum.tar.Z 

'* IRAF is distributed by the National Optical Astronomical Observatories, 
which are operated by AURA, Inc., under cooperative agreement with the 
National Science Foundation. 

ECLIPSE is a software package written by N. Devillard, which is avail- 
able at http://www.eso.org/projects/aot/eclipse 



For each science exposure in a given OB, we constructed a 
"sky" image to remove the sky and telescope background, and 
the dark+bias signal. We used a minimum of three and a max- 
imum of eight temporally adjacent frames depending on when 
the science exposure was taken during the OB sequence. In a 
"first pass," the individual sky frames were scaled to a com- 
mon median level equal to that of the science frame, averaged 
together with minimum-maximum rejection algorithm at each 
pixel, and subtracted from the science frame. Cosmic rays and 
hot pixels were identified using the "L. A. Cosmic" routine 
(van Dokkum 2001) and negative outliers at > 8<t below the 
mean residual background were added to the bad pixel lists. 
The sky-subtracted frames within each OB were aligned using 
integer-pixel shifts and coaveraged. From this intermediate 
combined OB image, an object mask was created by applying 
an appropriate threshold to identify astronomical sources. In 
a second "mask pass," the sky subtraction procedure was re- 
peated using the object masks to minimize the contribution of 
sources in the background computation. The mask-pass sky- 
subtracted frames were divided by the master flat field for the 
corresponding epoch of observation and band, and combined 
as before. Associated weight maps were generated, propor- 
tional to the total integration time of all frames contributing at 
each pixel. We normalized the final OB images to a 1 s inte- 
gration time and the weight maps to a maximum value of 1 . 

Non-linearity (raw counts exceeding 10000 ADU, for 
which ESO reports 99% linearity) is not a concern except for 
the very few brightest sources (< 3 per ISAAC field) in the 
frames obtained under the best seeing and transparency con- 
ditions. We did not apply any correction for non-linearity but 
excluded those stars with raw counts well above 10000 ADU 
in the flux calibration and PSF determination described below. 

In coaveraging the reduced frames, bad pixels and cosmic 
rays were masked out and no clipping or rejection algorithm 
was applied. We verified that no strongly deviant pixel, either 
persistent or as single event, was missed by comparing the OB 
images to versions obtained by median-combining the regis- 
tered sky-subtracted frames without masking and from the be- 
haviour of individual pixels throughout the stack of unaligned 
unmasked sky-subtracted frames. In several OBs, we noticed 
satellite tracks, moving targets (presumably Kuiper Belt ob- 
jects), and spurious features produced by reflection effects of 
bright sources falling on the edges of the detector array. We 
removed them with customized masking. We inspected the 
individual reduced frames to identify those with significant 
residual artefacts due to poor sky subtraction or instrumen- 
tal features. We implemented additional mask-pass routines 
in DIMSUM to improve the quality of the affected data, de- 
scribed in what follows. 

Bias residuals — Particular processing was required to ac- 
count for the time and illumination-dependent bias level of 
the detector The median-scaling applied in constructing the 
sky images accounts for the global intensity of the dark+bias 
but not for its spatial structure. Because the bias variations 
are strong along columns but negligible along rows, we re- 
moved the residuals by subtracting the median row by row 
in all sky-subtracted frames. In many OBs, the bias var- 
ied so strongly and rapidly that it led to spurious "numeri- 
cal ghosts" with shapes reminiscent of the masked sources. 
These features were most conspicuous in regions of steepest 
spatial gradients in the bias (lower parts of the bottom and top 
half of the frames) and during periods of steepest gradients 
in time (related to the illumination history). The average of 
the scaled sky frames did not represent accurately the local 
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background in the science frame and the effect was worsened 
by the object masking which reduced locally the number of 
frames used to estimate the background. The computed back- 
ground in the regions with objects masked in one or more of 
the sky frames thus differed significantly from that of adja- 
cent regions masked in none of the sky frames, imprinting 
positive ghosts for a locally underestimated background or 
vice-versa. We eliminated the numerical ghosts while con- 
structing the sky frames: in addition to the median-scaling, 
we accounted for an additive component representing the bias 
offsets between a science frame and each of its corresponding 
scaled sky frames by fitting a median to each row in the dif- 
ference maps. We applied this procedure to the 28 OBs where 
ghosts were clearly recognized by eye inspection, mostly in 
the Js and H bands. 

Residual large-scale gradients — For about half of the 
data, significant large-scale gradients remained in the sky- 
subtracted frames. They were caused primarily by highly 
variable sky and thermal background emission, and stray 
light. The fractions of OBs affected in each band reflect the 
background behaviour seen in Figure 92% in H, 46% m 
Ks, and 27% in 7s- We first attempted to remove these resid- 
uals by fitting two-piece cubic splines directly to the object- 
masked frames, along rows or columns depending on the gra- 
dient pattern. Owing to crowding in the MS 1054-03 field, 
the fits were however generally poorly constrained. We thus 
created template "flat images" by registering and coaverag- 
ing all sky-subtracted frames unaffected by residual gradients, 
taken from all OBs of the appropriate field and band to in- 
crease the S/N ratio. We subtracted the flat templates from 
the affected frames. This largely eliminated sources, leaving 
residual cores due to small mismatches in alignment and see- 
ing. We used these difference maps to create new masks and 
fit the two-piece cubic splines. Because the masked area was 
substantially reduced, the fits were much better constrained, 
allowing successful removal of the residual gradients. We 
performed this second-order sky subtraction for the 29 OBs 
where the gradients were strongest and apparent in the com- 
bined OB images. 

Instrumental features and other artifacts — A high spa- 
tial frequency pattern known as the "odd-even column effect" 
appeared in part of the data taken after the 2001 March re- 
aluminization (19 OBs). It is characterized by an offset in 
level between adjacent columns, occasionally further mod- 
ulated over a width of about 20 pixels and generally differ- 
ent between detector quadrants. When present, we corrected 
for this artifact by subtracting the median of the background 
along columns and for each quadrant separately in the object- 
masked sky-subtracted frames. 

Finally, high to moderate spatial frequency ripples presum- 
ably from moonlight reflections were present in the back- 
ground of 41 OBs, with important frame-to-frame variations 
in intensity, structure, and orientation. A semi-circular high- 
frequency time-dependent interference pattern was also obvi- 
ous in 3 OBs. The complexity and variability of these artifacts 
make them very difficult to correct for Since they generally 
averaged out in the combined OB images, we did not treat 
them in any particular way except by excluding the most af- 
fect frames which left noticeable patterns in the OB images. 

3.1.3. Rectification 

We corrected the OB images for the documented ISAAC 
geometrical distortion while simultaneously registering and 
interpolating the ISAAC data to the WFPC2 maps at 



0."1 pixer^ This "rectification" was done using the same 
algorith m as for the HDF-S data described by Labbe et aQ 
J2003ah . Briefly, the transformation corrects for distortion, 
aligns the ISAAC OB images (allowing for sub-pixel shifts), 
and applies the best-fit linear transformation determined from 
the positions of stars and point-like sources across the 7s -band 
images and the WFPC2 /g^-band mosaic. The images are re- 
sampled using third-order polynomial interpolation, and only 
once to minimize the effects on the noise properties. We 
rescaled the rectified OB images to ensure flux conservation. 
Since the ISAAC field distortions have small amplitude and 
thus negligible effects on the photometry, we used a constant 
flux scaling factor for the entire images. 

The rectification of the OB images instead of the individ- 
ual frames affects the point-spread function (PSF) because of 
the dithering and integer-pixel registration within each OB. 
Figure |3 plots the PSF FWHM measured in the rectified 
OB images versus the median of the FWHMs in the single 
frames of the corresponding OB. The effect is negligible at 
FWHM > 0."5 and increases towards lower FWHM up to a 
maximum of 25% at 0."3, the limit of undersampling. 

3.1.4. Flux Calibration 

We determined the instrumental photometric zero points 
(ZPs) using observations of standard stars. All stars were se- 
lected from the LCO/Palomar NICMOS Hst ( Persson et aQ 
1998) except for one taken fro m the UKJRT Faint Standards 
list (iCasali & Hawardeill 19921 observed on 1999 December 
18). We measured all fluxes and magnitudes for photometric 
calibration in a 6 "-diameter circular aperture, which proved to 
be adequate by monitoring the growth curve of the standard 
stars and of reference stars in the MS 1054-03 field. 

The standard stars were observed in sequences of five 
dithered exposures. We integrated the count rates of the stars 
within the 6"-diameter aperture in each frame of a sequence 
after sky subtraction and flat-fielding (using the appropriate 
flat field for the epoch of observation and band). We com- 
puted single-star fluxes by averaging, with 2(T-clipping, the 
count rates for each sequence. Typical uncertainties, taken 
as the median of the standard deviations of single-star fluxes, 
were of 1.7% (average of 2.1%). Using the calibrated mag- 
nitudes, we derived single-star ZPs. We identified as non- 
photometric the nights where single-star ZPs varied by more 
than twice the typical uncertainties (i.e. 3.4%), and excluded 
them from the subsequent analysis. For each of the remaining 
nights, we derived nightly ZPs by averaging the single-star 
ZPs. Because the nightly ZPs exhibit a large systematic in- 
crease after realuminization events, we computed period ZPs 
by averaging with 2(T-clipping the nightly ZPs for each of PI 
to P3. Nights rejected by the 2o'-clipping were added to the 
initial list of non-photometric nights. In total, 17 out of 36 
nights were photometric according to our criteria. 

The median night-to-night dispersion of the ZPs is 
0.024 mag and varies from 0.007 and 0.060 mag depending 
on the band and period. The average ZPs per period in 7s, H, 
and Ks increased by 0.340, 0.415, and 0.441 mag from PI to 
P2, and by 0.276, 0.255, and 0.188 mag from P2 to P3. We 
adopted as reference ZPs for our MS 1054-03 data those of 
P3, which are reported in Table|2l 

We did not account for atmospheric extinction or for colour 
terms. The calibration stars observations cover airmasses 
1 - 1 .6, at which 95% of the MS 1054-03 data were obtained. 
As for the MS 1054-03 stars count rates (Figure|3|l, the stan- 
dard stars data do not show any measurable correlation with 
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airmass. The ISAAC H and Ks filters m atch well those used to 
establish the "LCO" standard system of lPersson et al. I (Il998h 
while the 7s filter has a slightly redder effective wavelength 
than the LCO J filter For a colour term ^{J-K)i^co, the 
colour coefficient 7 is theoretically e xpected to be < 0.01 mag 
in H and K^, and -0.04 mag in 7s (lAmico et al. 2002) but 
these values have not been verified experimentally yet. The 
ZPs derived for the individual stars do not reveal any depen- 
dence on J-K colour over the range 0.24-0.46 mag covered 
by the standard stars set, with a scatter much larger than the 
expected colour effects. 

We calibrated the MS 1054-03 rectified OB images using 
a set of 4-5 bright, isolated, and unsaturated stars for each 
field. We first considered only the OBs taken during photo- 
metric nights. We integrated the stars count rates in a 6"- 
diameter aperture and applied the nightly ZPs to obtain indi- 
vidual magnitudes. We adopted as calibrated magnitudes the 
2cr-clipped averages of all individual measurements. We then 
included all OBs and computed for each star and OB the ap- 
propriate scaling factor based on the star's magnitude and the 
adopted ZPs of Table |2] We calibrated the OB images indi- 
vidually using the 2cr-clipped average of the scaling factors of 
all stars per OB. 

For each reference star in the MS 1054-03 field, we com- 
puted the 1(7 dispersion of their calibrated magnitudes in the 
photometric OBs. The median of the dispersions of all ref- 
erence stars is 0.030, 0.032, and 0.025 mag in 7s, H, and K^, 
respectively. We adopt 0.03 mag as the internal accuracy of 
the photometric calibration of the NIR MS 1054-03 data set. 

3.1.5. Fina I Mosa icking 

We produced the final 7s-, H-, and /Ts-band images of 
MS 1054-03 in two steps: we first combined all OB images 
per field and then combined the four field maps into the large 
mosaic. We proceeded in this way because of systematic dif- 
ferences in PSF of the OB images between the data sets of 
each field, which would be problematic for the overlap areas. 

To create the separate field maps, we weight-averaged the 
calibrated rectified OB images per band, without clipping or 
rejection algorithm. The weight applied at each pixel vvpix 
was proportional to the integration time fint.pix, scaled by the 
inverse of the pixel-to-pixel variance of the background noise 
in each image CTqb ™d of the squared FWHMob of the OB 
PSF: 

(1) 



Wnix OC 



^iiit.pix 



ao^BFWHMo^B 

The fint.pix term corresponds to the weight maps of the OBs 
multiplied by their respective total integration time, and thus 
accounts for pixels that were excluded because bad/hot or af- 
fected by other artefacts. The background noise was estimated 
on source-free areas of the calibrated but pre-rectified OB im- 
ages to avoid the spatially-correlated noise introduced by rec- 
tification. The CTqb FWHMqb term represents the noise within 
an aperture of the size of the seeing disk (assuming purely 
Gaussian noise) so that the weighting scheme optimizes the 
S/N ratio for point sources. We created weight maps for each 
field by averaging those of the individual OBs applying the 
same weighting terms (i.e., given by Eq. [0. These maps, 
hereafter denoted wp, give the weighted exposure time per 
pixel. They reproduce accurately the spatial variations of the 
background noise within each field, as indicated by the con- 
stant pixel-to-pixel rms over the "noise-equalized" maps ob- 
tained by multiplying the images with the squared-root of the 



associated weight maps. 

The sky subtraction algorithm in DIMSUM and our ad- 
ditional optimization routines have introduced small nega- 
tive biases in the background visible as orthogonal stripes at 
P. A. w 30°, and as dark areas around large sources and in 
crowded regions. They are due to an overestimate of the back- 
ground levels because of insufficient source masking. The ob- 
ject masks were defined from the shallower OB images, such 
that emission in the extended PSF wings around sources was 
not fully accounted for and very faint sources were missed. 
Since this effect is systematic in all OBs, it is enhanced and 
becomes apparent in the deeper field maps. To reduce these 
residuals, we rotated a copy of the field maps, made new ob- 
ject masks using SExtractor (as described in § |5}, fitted the 
background in the unmasked regions with bi-cubic splines 
along rows and columns, rotated the fits back, and subtracted 
them from the field maps. 

Our aim was to produce final mosaics as uniform as possi- 
ble in effective angular resolution and optimize the S/N ratio 
in the overlap areas. We determined the effective PSF of the 
field maps by averaging sub-images centered on bright, iso- 
lated, unsaturated stars (normalized to the stars' peak flux). 
Table [T] lists the FWHM of the best-fit Moffat profile to the 
PSFs. For each band separately, we convolved the three high- 
est resolution fields to bring all fields to a common resolution. 
We computed the kernels with a Lucy-Richardson deconvolu- 
tion algorithm. We quantified the accuracy of the PSF match- 
ing from the curve-of-growth analysis of the convolved PSFs, 
measured from the same stars in the smoothed field maps. For 
each band, the fraction of enclosed flux agrees to better than 
1% in circular apertures with diameters in the range l"-2", 
and within 4% for diameters between 0."7 (w 1.5 x FWHM 
of the PSF of the smoothed field maps) and 6" (the reference 
aperture for photometric calibration). 

To optimize the S/N ratio in the overlap areas, we applied 
a similar weighting scheme as in Eq. ^ with however some 
differences. We omitted the FWHM-dependent term since it 
is identical among the PSF-matched field maps. We needed 
to account for the effective background noise levels properly. 
Using the pixel-to-pixel background rms determined over the 
whole field images resulted in significant global differences 
between the four fields in the noise-equalized mosaics. This 
indicates that the relative rms-weighted total integration times 
between the fields does not reflect the differences in effec- 
ti ve n oise. We attribute this to the systematics discussed in 
§ 13.11 which affect in different proportions and with varying 
severity the data of each field. Thi s is evident in the results 
of the noise analysis described in § 14.41 where the deviations 
from the behaviour for (uncorrected) pure Gaussian noise of 
the real background rms as a function of aperture size is no- 
ticeably different between the fields. For the inverse variance 
term of the weights, we used the effective noise derived from 
this analysis for an aperture diameter of 1.5 x FWHM. This 
led to the most uniform noise-equalized mosaics in terms of 
background pixel-to-pixel rms, and again optimizes for point 
sources. 

We created three sets of weight maps for the mosaics, all 
based on the fields' weight maps wp. We made one set by 
averaging the wp maps, applying the same weights as for 
the images. This set accounts for the variations of effective 
background noise over the entire mosaic relevant for photom- 
etry and source detection based on a S/N threshold ("effective 
weight maps," with values denoted by Weff). We made a sec- 
ond set weighting only by the cumulative number of frames 
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contributing to the flux (i.e., with only the fint,pix and no vari- 
ance term), as record of the relative total integration time per 
pixel ("exposure time weight maps," Wexp). For the third set, 
we simply averaged the wp maps assigning equal weights so 
that in the mosaic, the maximum weight in the deepest part of 
all four fields is identical. This set is appropriate for detect- 
i ng so urces using a surface brightness threshold as we do in 
§ l5.1l ("adiusted weight maps," Wadj). We renormalized all the 
weight maps for the mosaics to a maximum of unity. 

3.2. FORSl Data Reduction 

3.2.1. Bias, Flat Field, Sky Subtraction, and Cosmic Rays 
Rejection 

The bias subtraction and flat-fielding of the optical FORS 1 
imaging were carried out using pipeline reduction tasks de- 
veloped fortheFORSinstruments within the ESO-MIDAS 
system jWarmelslfT99l . Bias subtraction was performed us- 
ing the nightly master bias frame scaled to the average level 
of the science frame in the pre- and over-scan regions. Flat- 
fielding was performed with the nightly master sky flat image 
obtained from series of exposures taken during twilight. 

We further reduced the science frames within IRAF as fol- 
lows. We subtracted a constant sky level corresponding to the 
mode of each frame. In one f/-band OB, there was a residual 
offset in background level between the lower and upper half 
of the array, which we eliminated by adding to the top half the 
difference in median over 20 rows below and above the mid- 
array row 1024. We identified cosmic rays with L.A.Cosmic 
ivan Dokkuml200lh and negative bad pixels by applying a 3 cr 
threshold below zero. We updated the resulting masks to flag 
out additional artefacts identified by visual inspection (e.g., 
occasional bad rows and columns, transient positive source- 
like features and broad faint tracks, edges of the low sensitiv- 
ity "holes" intrinsic to the FORS 1 array not adequately flat- 
fielded out). 

We removed residual large-scale patterns (at a few percents 
of the subtracted constant sky level) using a background map 
gener ated with the SExtractor program ( Bertin & ArnoutsI 
|l29^ run on the initial sky-subtracted frames. We specified a 
grid mesh size of 64 pixels and median filter size of 10 meshes 
for the background estimation, which produced the best re- 
sults. It is unclear whether these residuals were due to flat 
field inaccuracies or spatial structure in the background. Ei- 
ther division by or subtraction of the SExtractor background 
map reduced the residuals to < 1%. The two methods led 
to essentially indistinguishable results, with differences in re- 
maining gradients well below 1 % and in photometry of sev- 
eral stars in the field smaller than 0.002 mag. As reported on 
the ESO FORS 1+2 Web pages, flat fielding with the master 
sky flat can leave large-scale gradients of a few percents, so 
we adopted the sky-subtracted frames divided by the SExtrac- 
tor background maps. 

3.2.2. Flux Calibration 

We determined the photometric ZPs from observations of 
standard stars taken during the same nights as the MS 1054- 
03 imaging. These were available for the single OB in the V 
band, and two OBs in each the B and U bands. The stars 
were selected from the catalogue of iLandoltl ( 119921) . Sev- 
eral exposures were obtained through each night, and each 

Details of the pipeline reduction for FORS 1+2 can be found at 
Pittp ://www.eso.org/observing/dfo/quality/FORS/pipeline/pipe_reduc.html| 



exposure includes multiple standard stars. The stars data 
were bias-subtracted and flat-fielded in the same way as the 
MS 1054-03 data. We determined single-star ZPs from in- 
dividual star counts, and accounted for the non-negligible 
effects of airmass and colour terms using the atmospheric 
extinction and colour coefficients for the corresponding pe- 
riod of observations reported on the ESO/FORSl Web pages. 
We computed the nightly ZPs by averaging with 2cr-clipping 
the single-star ZPs. The \a dispersions range from 0.013 to 
0.049 mag, depending on the night and band. We also com- 
puted single-frame ZPs as the 2(T-clipped average of the mul- 
tiple stars ZPs per frame. Frame-to-frame variations in ZPs 
are smaller than the single-star ZPs dispersion and there is 
no trend with time, indicating that the nights were photomet- 
ric. As for the ISAAC data, we used a circular aperture of 
diameter 6" for the photometric calibration. We verified that 
no measurable systematic variations with airmass or colours 
were left. The adopted ZPs as well as the extinction and 
colour coefficients used are given in Table |2 

For the MS 1054-03 data collected on nights without ap- 
propriate standard stars observations (one OB in B and three 
OBs in U), we proceeded as for the ISAAC data taken during 
non-photometric nights. We first determined the magnitudes 
of five bright, isolated, unsaturated stars from the subset of 
calibrated MS 1054-03 frames, averaging with la clipping 
the individual measurements (the mean of the Icr dispersions 
are 0.013 mag in B and 0.036 mag in U). We then computed 
for each star and frame the appropriate scaling factor based 
on the stars' magnitudes and adopted ZPs. We calibrated the 
frames individually using the 2o'-clipped average of the scal- 
ing factors of all stars per frame. We estimate the internal 
accuracy of the photometric calibration of the FORS 1 optical 
MS 1054-03 data set to be 0.015 mag in V, 0.02 mag in B, 
and 0.04 mag in U . 

3.2.3. Combination and Rectification 

We aligned all calibrated frames allowing for integer-pixel 
shifts. We then co-averaged all registered calibrated frames 
for each band, applying the weighting scheme of Eq.^with 
the various terms determined from the individual frames. This 
was done before rectification to the WFPC2 maps. Because 
of the smaller number of frames and the poorer seeing con- 
ditions for the FORS 1 observations (compared to the ISAAC 
observations), the smearing due to the combination of dither- 
ing, uncorrected field distortion, and sub-pixel misalignment 
are not critical. The co-averaging, rectification, and creation 
of the weight maps were done in the exact same way as for 
the ISAAC field maps. 

4. FINAL IMAGES 
4.1. General Properties 

The final reduced images have a pixel size of Q'.'l and are 
aligned relative to the WFPC2 I%n mosaic, with North ori- 
ented near ly up at 0.°78 anticlockwise from the vertical axis 
(see § 14. 31 for the exact astrometry). The images are all nor- 
malized to instrumental counts per second. The fully reduced 
images as well as the raw data are available electronically on 
the FIRES Web site'^. 

The ISAAC NIR mosaics and FORS 1 optical maps are shal- 
lower towards the outer edges, reflecting lower total exposure 
times because of the dithering process in the observations. 

' ^ http ://w w w. strw. leidenuniv. nl/~fires 
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In the NIR mosaics, the largest weights are reached in the 
regions where the adjacent fields overlap (12% of the total 
area mapped with ISAAC). In addition, the combination of 
different integration times and noise levels results in signif- 
icant global offsets in weight between the four fields. The 
largest effective weights in the central parts of the fields are 
Weff « 0.9 for all NIR bands (in Fl for J, and K^, and F4 for 
H) compared to the normalized maximum of 1 reached in the 
overlap regions. Field F2 has the lowest maximum effective 
weight, reaching Weff = 0.64, 0.5, and 0.39 in J^, H, and K^, re- 
spectively. This can bias area selection when applying weight 
criteria. To take advantage of the field size, it may be more 
appropriate to use the adjusted weight maps where the max- 
imum weights are equal in all fields (see § 13. U or to treat 
the fields individually using their respective weight maps wp. 
The area in the /Q-band mosaic with Wadj > 0.8, 0.25, and 
0.01 (corresponding approximately to wp > 0.95, 0.30, and 
0.01 in at least one of the four ISAAC fields) is 16.3, 26.2, 
and 28.7 arcmin^. As quality cut for photometry, we adopt 
the 23.6 arcmin^ area defined by Wadj > 0.25 (approximately 
wp > 0.3 in at least one ISAAC field) in each NIR band, and 
w > 0.3 in each optical band. 

The noise-equalized /Q-band mosaic (i.e. the /Q-band map 
multiplied by the squared root of its Weft weight map) is shown 
in Figure [2 It illustrates the richness in faint sources and 
the spatial uniformity in normalized background rms achieved 
with the effective weight map. Figure |5] shows the RGB 
colour composite image of the central 4' x 4' of the field, con- 
structed from the /814, J^, and maps. The three maps have 
been convolved to an angular resolution of FWHM = 0."69 
to match the broadest PSF of all images (in the U band; see 
§ 14.21 below). The colour image is displayed with a linear 
stretch adjusted to enhance the fainter objects. The variety of 
colours is clearly visible over a wide range of brightnesses, 
even at the faintest levels. This indicates that the maps reach 
comparable depths, and reflects the combination of different 
rest-frame colours and different redshifts of the sources. 

4.2. Image Quality and PSF-Matching 

We characterized the PSF in the final combined maps us- 
ing bright, isolated, and unsaturated stars (the sets were partly 
different in the various bands and included eight to thirteen 
stars). The PSF in the ISAAC and FORSl images is stable 
and symmetric over the entire field of view, with a Gaussian 
core profile and wings consistent with a Lorentzian profile. 
We adopted the average of the peak-normalized star profiles 
as the final PSF in each band. From Moffat profile fits, the 
PSF FWHM in the J,-, H-, and K,-hand mosaics is 0"48, 
0"46, and 0."52, with ellipticities of 0.02, 0.04, and 0.03, re- 
spectively. For the optical bands, the PSF FWHMs and ellip- 
ticities ai-e 0."69 and 0.06 in U, 0"57 and 0.04 in B, 0"65 and 
0.02 in y, 0."21 and 0.02 in Veoe, and 0."22 and 0.02 in ki^. 
The FWHM of the individual stars are within w 10% of the 
final PSF FWHM in each of the ISAAC and FORSl images, 
and within w 30% for the WFPC2 images. The PSF FWHMs 
are reported in Tabled 

For consistent photometry across all bands, we convolved 
all images to a common angular resolution corresponding to 
that of the f/-band map, which has the poorest resolution 
(FWHM = 0"69). The convolution kernels between the U- 
band PSF and that of the other bands were computed using a 
Lucy-Richardson deconvolution algorithm. From the curve- 
of-growth analysis of the convolved PSFs, the fraction of en- 
closed flux agrees to < 2% in circular apertures with diame- 



ters between 1" and 2", relevant for the colour measurements 
(see § 15. 2> . The accuracy of the PSF-matching is < 4% for 
diameters in the range l"-6". 

4.3. Astrometry 

Accurate registration among the various maps is important 
for correct cross-identification of sources, precise colour mea- 
surements, and proper comparisons of the morphology in dif- 
ferent wavebands. We verified the relative alignment using 
20 bright stars and compact sources throughout the field. The 
median rms variation in position of individual objects across 
the bands is 0."07 (or 0.7 pixel in the final images), with a 
range from 0"025 to 0"20. Relative to the positions in the 
WFPC2 /814 or Veoe maps, there is no systematic trend in the 
position offsets in any of the bands. The amplitude of the 
offsets in the J^, H, and Ks mosaics has a Ict dispersion of 
w 0"07, corresponding to a fraction 0.47 of a pixel at the 
original pre-rectified scale of the ISAAC data (0"148 pixel"'). 
For the U, B, and V maps, the offsets dispersion is w 0."12, or 
0.6 pixel at the original FORSl scale (0"20). 

We verified the absolute astrom etry using s tars from the 
USNO ca talogue (v. A2. 0, Monet et al.ll1998h and UCAC2 
catalogue fZacharias et al. 2004). We solved for the positions 
determined from the WFPC2 /814 mosaic, our reference image 
for the relative registration. The astrometric solution is 

A R. A. =-xs cos(ip) -ys sm(ip) + 208 ."566 
ADEC = -xi sin(v3) + yicos((^)- 194."435 (2) 

where x and y are the pixel coordinates in the final images, 
for which the precise spatial scale is s = 0."0996 pixel"'. The 
orientation is such that North and East coincide nearly with 
the vertical and horizontal axes, respectively, with a small ro- 
tation angle of ip = 0.°77754 counterclockwise (following the 
usual astronomical convention. East is 90° counterclockwise 
from North). The offsets AR.A. and A DEC are expressed 
in arcseconds and relative to the reference coordinates Q!2oo{): 
10''57"00?07, (520()o: -03°37'36."19. We estimate that the ab- 
solute astrometry is accurate to < 0"5. 

4.4. Background and Limiting Depths 

In the raw data, the background noise properties are ex- 
pected to be well described by the dispersion of the signal 
measured in each pixel since both the Poisson and read-out 
noise should be uncorrected. For such uncorrected Gaus- 
sian noise, the effective background rms for an aperture of 
area A is simply the pixel-to-pixel rms a scaled by the linear 
size = VA of the aperture, a(N) = Na. Instrumental fea- 
tures, the data reduction, and the rectification, combination, 
and PSF-matching procedures have added significant system- 
atics and correlated noise in our final data. The determination 
of the limiting depths and photometric uncertainties depends 
critically on the accurate description of the noise properties of 
the images. For this purpose, we followed the same e mpirical 
approach as for the HDF-S data iLabbe et al. 12003' a^. 

We derived the function a{N) directly from the data by mea- 
suring the fluxes of 2800 non-overlapping circular apertures. 
We placed the apertures at random in the maps but so as to 
avoid all pixels associated with sources detected in the Ks 
band and within their ij,(Ks.ab) = 24.8 magarcsec"^ isophote 
(the surface brightness threshold applied for source detection; 
see § |5}. We used the same positions for all maps, and re- 
peated the measurements with different aperture diameters 
ranging from 0."7 to 3". Because the noise properties vary 
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importantly between the four fields of the NIR mosaics, we 
analyzed them separately (each field has about 700 aperture 
positions). We applied the procedure on the PSF-matched im- 
ages used for the photometry, allowing us to compute realistic 
photometric uncertainties. We also ran the procedure on the 
set of non-concolved WFPC2 and ISAAC images to estimate 
more representative intrinsic depths, which we can also better 
compare with those of the HDF-S data (see below). 

For each aperture size, the distribution of fluxes is symmet- 
ric around the background value (about 0), and is well repro- 
duced with a Gaussian profile. We thus determined the back- 
ground rms variations from the FWHM of the best-fit Gaus- 
sian to the observed distribution for each aperture size. While 
at any fixed spatial scale, the noise properties over the im- 
ages are consistent with a pure Gaussian, the variations with 

deviate appreciably from Gaussian scaling, being system- 
atically larger and increas ing non-linearly. As for the HDF- 
S data jLabbe et~ 1 12003 a). the real noise behaviour in the 
MS 1054-03 data set can be well approximated by the poly- 
nomial model: 

Naiiai + biN) 
ai(N) = ^ , (3) 

where / designates the band (and the field for the ISAAC data) 
and the weight term accounts for the spatial variations in noise 
level related to the exposure time. The coefficient a represents 
the effects of correlated noise, which are dominated by the 
resampling of the images to the WFPC2 pixel scale and, for 
the PSF-matched data, the convolution applied. The curvature 
in the cr,(A^) relations quantified by the coefficient b indicates 
a noise contribution that becomes increasingly important on 
larger scales. It presumably originates from objects below the 
detection level, instrumental features, and residuals from the 
sky subtraction and flat-fielding. 

The noise behaviour is qualitatively the same in all of the 
images. Figure |6l shows the results for the four fields of the 
PSF-matched /Tj-band mosaic, with the background rms mea- 
surements for the various apertures and the best-fit polyno- 
mial models (Eq.|3} compared to the expected linear relation- 
ship for uncorrelated Gaussian noise. The flux distributions 
from which the background rms is derived are also plotted for 
selected aperture sizes. With the empty areas defined from the 
/Cj-band map, sources that are undetected in but detectable 
in the other bands will contribute to the aperture flux distri- 
butions. This is reflected in the distributions for bands other 
than Ks, more obviously for the optical bands, with a small 
excess of positive flux measurements. However, the bias in- 
troduced in the results is very small because the positive tail 
affects very little the Gaussian fits. 

Table |3l summarizes the results for the PSF-matched im- 
ages, along with the corresponding la background-noise lim- 
iting magnitudes. Table |4]gives the limiting magnitudes for 
the rectified but unconvolved images. The limiting mag- 
nitudes refer to a "point-source aperture" with a diameter 
d = 1.5 X FWHM that maximizes the S/N of photometric 
measurements in unweighted circular apertures of point-like 
sources. Compared to the HDF-S data set, the MS 1054-03 
field limiting magnitudes are 0.5- 1.9 mag brighter, depend- 
ing on the band. This comparison is rather crude because it 
does not account for the partly different bandpasses in the op- 
tical (U and B in particular) or for the different angular reso- 
lution between the data sets. 

For the ISAAC NIR data, the J,-, H-, and K.-hmd PSF- 
matched maps are w 0.8 mag shallower than the HDF-S im- 



ages. The latter were matched to the resolution of the ISAAC 
//-band map of F'WHM = 0"48 and the point-source aper- 
ture has d = '.'1 ( iLabbe et al. 112003 a'). The intrinsic resolu- 
tion of the MS 1054-03 NIR rectified but unconvolved mo- 
saics is similar, enabling a better comparison. The limiting 
magnitudes derived from this set are w 0.65 mag brighter 
than for the HDF-S. The differences expected based only on 
the integration time are of « 0.9 mag for similar observing 
conditions (seeing, atmospheric transmission, and night-sky 
OH and continuum emission). The depths achieved for the 
MS 1054-03 field exceed slightly these expectations. Vari- 
ous factors can explain this performance. The most impor- 
tant one is probably that about half of the MS 1054-03 data 
were taken after the 2001 March realuminization (period P3), 
when the effective sensitivity had increased by w 20% - 80% 
compared to earlier epochs since the end of 1999 (PI and P2) 
during which the HDF-S obs ervations were carried out (§ 12.21 
see also ILabbe et"aL ,l2003ah . 

We also analyzed the distributions of flux ratios between 
bands. The results obtained for ratios measured in random 
pairs of apertures (i.e., not necessarily spatially coincident) 
correspond closely to the predictions based on Eq. |3l assum- 
ing the noise is uncorrelated between bands. In contrast, 
the scaling relations derived for ratios within registered aper- 
tures for several band pairs lie systematically below such pre- 
dictions on scales ofA^> 10-15 pixels (although still well 
above expectations for pure uncorrelated Gaussian noise). 
This suggests that the background noise pattern is correlated 
between the bands, as was also seen in the HDF-S data set 
("Labbe et al. "2003a'). The effects are most important and up 
to « 20% among the maps obtained with the same instrument. 
This is likely due to the large-scale residuals remaining after 
the data processing (e.g., from the bias, sky, and flat-field), 
which are similar for a given instrument. Interestingly, this 
effect is also present at the 10% level for ratios involving the 
FORSl B or y bands and the WFPC2 Va)6 or bands, as 
well as between the WFPC2 /814 and ISAAC bands. This is 
probably the signature of sources undetected in /Tj but clearly 
present in the other bands contaminating the "empty" back- 
ground flux measurements. Confusion noise was argued in 
the case of the HDF-S data but is unlikely for the shallower 
MS 1054-03 maps. The implications of this analysis is that 
our uncertainties on colours are probably slightly overesti- 
mated because we propagated those of each bandpass neglect- 
ing the cross-correlation terms. We did not attempt to account 
for these effects as they are much smaller than the difference 
between uncorrelated Gaussian scaling and our empirical re- 
lationships. 

5. SOURCE DETECTION AND PHOTOMETRY 

Our aim was to construct a /Ts-band selected source cat- 
alogue as uniform as possible over the entire MS 1054- 03 
field, optimized for point-like sources, and with reliable de- 
tections and optical-NIR photometry. Because of the differ- 
ent noise properties between the four ISAAC fields, detec- 
tion with a constant S/N threshold would result in different 
limiting magnitudes and limiting surface brightnesses from 
one field to the other We therefore applied a constant surface 
brightness threshold fiiKs) instead. We compromised between 
completeness and reliability by choosing a ii(Ks) criterion so 
as to keep the contamination by spurious (faint) sources be- 
low ^ 5% in all fields. The general procedure for source 
detection and photo metry is analoguous to that followed by 
ILabbe et al. I ll2()()3al) in making the HDF-S catalogue. 
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5.1. Detection 

We carried out the s ource detection using th e SExtractor 
software version 2.2.2 (Bertin & Arnouts 1996). The objects 
were detected in the rectified but non PSF-matched A'j-band 
map. To optimize for point-like sources, we filtered the de- 
tection map with a Gaussian kernel of FWHM = 0."52 ap- 
proximating well the core of the A's-band PSF. The detec- 
tion criterion was that at least one pixel be brighter than 
M(^s,ab) = 24.8 mag arcsec"-. This corresponds to w Scr of 
the background noise in the filtered /Q-band mosaic over the 
two fields with intermediate depth (F3 and F4). The threshold 
adopted is rather conservative but ensures that the spurious 
source contamination is small. 

The resulting catalogue contains 1858 objects. We esti- 
mated the fraction of false detections by running SExtractor 
with the exact same parameters on the inverse detection map; 
it is < 1% over the entire mosaic and < 3% in the shallowest 
field, F2. SExtractor identifies and separates blended sources 
based on the distribution of the /Tj-band light in the detec- 
tion map; 18% of the sources were flagged as blended. We 
specified deblending parameters so as to avoid oversplitting 
of objects with relatively complex morphologies. 

The spatial filtering affects the detection process and the 
isophotal parameters. Our filtering with a kernel reproduc- 
ing the /Q-band PSF introduces a small bias agai nst faint ex- 
tended sources. As for the HDF-S catalogue (see Lab be et al. I 
l2003ah . we preferred not to combine multiple catalogues ob- 
tained with different kernel sizes because their subsequent 
merging is complicated and subjective, and does not lead to 
a substantial gain in sensitivity for larger objects. The PSF- 
like kernel we adopted is fully adequate for our purposes since 
most of the faint sources detected are compact or unresolved 
in the NIR and it reduces the problems due to blending when 
convolving with larger kernels. 

5.2. Photometry 

We used SExtractor in dual-image mode, which allows for 
spatially accurate and consistent photometry across all bands. 
While the sources were extracted and the isophotal parame- 
ters determined from the filtered /Q-band map, the photometry 
was performed on the PSF-matched maps and within the same 
set of apertures in all bands for each object. Of the 1858 K^- 
band detected sources, 1663 have full eight-band photometry. 
We used circular apertures with 30 different diameters ranging 
from d = 0."7 to 3". We also used isophotal apertures defined 
by the detection threshold of iiiKf^.AB) = 24.8 mag arcsec"^ 
and autoscaling apertures, which, following iKronI (il98 0). 
are elliptical apertures scaled based on the first moments 
of the filtered TiTs-band light profile (the "MAG_ISO" and 
"MAG_AUTO" computed by SExtractor). Adopting the no- 
tation of Labbe et al. (2003a), we will refer to these aper- 
tures as APER(d), APER(ISO), and APER(AUTO), respec- 
tively. We also defined customized "colour" and "total" aper- 
tures APER(COLOUR) and APER(TOTAL). These were in- 
troduced for the HDF-S catalogue to allow detailed control of 
the photometry relying on well-defined criteria and simple but 
robust treatment of blended sources. We considered a source 
as b lended if the SExtractor fl ags "blended" or "bias" were set 
(see lBertin & Arnoutsll996h . 

The colour aperture is optimized for accurate and consis- 
tent colours, and relies on the /Tj-band isophotal aperture of 
each object. The criterion to select the appropriate colour 
aperture uses the equivalent circularized isophotal diameter 



diso = 2 -\/Aiso/7r. For isolated sources, 

r APER(ISO) 1 ."0 < Jiso < 2."0, 

APER(COLOUR) = <^ APER(1"0) di,o < 1"0, (4) 

iAPER(2"0) £/iso>2"0. 

For blended sources, 

C APER(£/i,o/i) 1 "0 < £/iso/i < 2"0, 

APER(COLOUR) = <^ APER( 1 ."0) di.Js < 1 ."0, (5) 

iAPER(2."0) di,Js>2'.'0, 

where s is the factor by which we reduce the size of the 
circular apertures centered on blended objects to minimize 
flux contamination from the close neighbour(s). The opti- 
mal value of s needs be determined from experimentation, 
and we found that the same 5 = 1 .4 as for HDF-S was ap- 
propriate for the MS 1054-03 field data. The smallest colour 
aperture aUowed, APER(1"0), has a diameter 1.5 x FWHM 
of the PSF-matched images. This maximizes the S/N of flux 
measurements for point sources in unweighted circular aper- 
tures and avoids smaller apertures for which the photometry 
becomes less reliable. The largest colour aperture considered, 
APER(2"0), was selected to avoid the uncertainties of very 
large isophotal areas and to reduce possible contamination of 
the fluxes by adjacent sources undetected in the band but 
present in other bands. 

The total apertures were used to derive the integrated fluxes 
of the objects in the /Tj-band. Isolated and blended sources 
were again treated differently: 



/ APER(AUTO) isolated sources 



(6) 



/\riiK(^ 1 u i^L.) - <^ APER(COLOUR) blended sources.' 

Equivalent circularized apertures were also defined for the au- 
toscaling Kron-like apertures with t/auto = 2 \/A^mo/t^ and for 
the total apertures with (i,o, = 2 ^Atot/7r. To compute the total 
fluxes and magnitudes, we applied an aperture correction to 
the measurements in APER(TOTAL). The correction factor is 
the ratio of the fluxes enclosed within circular apertures of di- 
ameters 6" (used for the photometric calibration) and c/tot, in- 
ferred from the grow th cu rve of the /Q-band PSF constructed 
from bright stars (§ 14. 2> . This correction is significant for 
most objects, reaching 0.58 mag for faint point-like sources 
with flux measured in a c/ = l."0 aperture. 

Figure0compares various methods of estimating integrated 
As-band magnitudes. For point-like sources, the total magni- 
tude is expected to correspond to the flux in the point-source 
APER(1"0) scaled by the aperture correction. The isophotal 
and autoscaling apertures tend to decrease in size at fainter 
magnitudes, leading to more severe underestimate of the total 
magnitudes due to extended emission intrinsic to the sources 
or in the PSF wings outside of the apertures. For isopho- 
tal magnitudes, the turnover sets in at Kl°)^g w 22.5 mag, 

with the difference [/:s(l"0)-0.58] growing rapidly 

and reflecting mainly the shrinking APER(ISO) defined by 
the constant surface brightness threshold. Magnitudes within 
APER(AUTO) come closer but still miss significant flux in 
the PSF wings at 7i:'°|g > 24 mag. With the method based 
on Eq.|6l we best recover the integrated magnitudes expected 
for faint point-like sources. Our treatment of isolated and 
blended sources also reduces the dispersion of Kl°^ around 
^:s(l"0)-0.58. 

We derived the la photometric uncertainties from Eq. |3l 
using the a,- and bt coefficients given in Table|3] the linear size 

= VA of the aperture considered, and the average weight 
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within the aperture. These uncertainties properly account for 
the actual background noise in the maps influencing directly 
the flux within the various apertures. They may o vere stimate 
the uncertainties of colours by a small amount ("§ I4.4> . From 
the variations in photometric ZPs determined from standard 
stars or inferred from stars within the MS 1054-03 field of 
view (§ 13. H and § I3.2> . we estimate that uncertainties of the 
absolute flux calibration are « 0.05 mag. The photometry 
may further be affected by possible biases from the sky and 
background subtraction procedure, or surface brightness bi- 
ases; these additional sources of errors are not included in our 
uncertainties. 

6. REDSHIFTS 

For all objects in the catalogue, we derived photomet- 
ric redshifts (znh) using a t echniqu e described i n detai l 
by Rudnicketal. (2001, 200j and iLabbe etal.l (Hooli)- 
Briefly, the algorithm involves linear combinations of red- 
shifted spectral templates of galaxies of various types (from 
normal ellipticals and spirals to magellanic irregulars and 
young starbursts). The mean absorption by neutral hydrogen 
from the intergalactic medium along the line of sight is ac- 
counted for (based on Madau 1995). The Zph uncertainties are 
estimated from Monte-Carlo simulations, accounting for un- 
certainties in the fluxes as well as template mismatch. A key 
feature of the technique is that the derived uncertainties reflect 
not only the la confidence interval around the best solution 
but also the presence of secondary solutions with comparable 
likelihood. 

Spectroscopic redshifts (zsp) have been measur ed for over 
400 galaxies in the MS 1054-03 field jTranet al. 1 1999. 2003 
Ivan Dokkum et a l^ 2003. 2004. S. Wuyts et al., in prepara- 
tion), of which about 330 are included in our /Tj-band selected 
catalogue. The majority (w 85%) lie at Zsp < 1, mostly mem- 
bers of the MS 1054-03 cluster at z = 0.83. Of the remaining, 
32 galaxies have 1 < Zsp < 2 and 1 1 have Zsp > 2. The much 
larger number of Zsp determinations in the MS 1054-03 field 
compared to HDF-S allows a more robust assessment of the 
accuracy of our photometric redshift technique. 

Figure|8]compares the Zph and Zsp for galaxies with spectro- 
scopic determinations, and shows the Zph and Zsp distributions 
of sources at Kl°}^j^ < 24.8 mag, for which s)N(KI%) > 5. 
Overall, the accuracy of the photometric redshifts is dz/(l + 
Zsp) = (|zph -Zsp 1/(1 + Zsp)) = 0.078. In the redshift ranges 
z < 1, 1 < z < 2, and z > 2, dz/(l+Zsp) = 0.074, 0.127, and 
0.040, respectively. At 1 < z < 2, the Zph's are less well con- 
strained, as reflected by the larger Icr confidence intervals. 
This is mainly because the fitting procedure there relies heav- 
ily on the presence of the Balmer/4000 A break in the SEDs 
(characteristic of > 10** yr old stellar populations), and the 
break falls in the large gap between the /814 and bandpasses. 
At z > 2, the agreement between Zph and Zsp is best, although 
the number of galaxies with available Zsp in this redshift in- 
terval is still small. Our Zph's slightly underestimate the true 
redshifts by Zph-Zsp ~ 0.1 on average. In the rest of the pa- 
per, we will use the spectroscopic redshift whenever available, 
otherwise the photometric redshift. 

We identified stars using a combination of three criteria. 
The main criterion was that the raw for fits using stellar 
SEDs was lower than that of the best-fit combination of galaxy 
SEDs from the Zph determination. The stellar templates are 
synthetic SEDs of main-sequence stars with effective temper- 
atures from 3000 to 10 000 K computed with th e NEXTGEN 
model atmospheres of iHauschildt et aO il999l) . We did not 



consider cooler dwarfs or evolved stars (giants, supergiants) 
because non local thermodynamical equilibrium (non-LTE) 
effects and extensive molecular absorption become important 
and model atmospheres still suffer from large uncertainties in 
these regimes. Moreover, adequate empirical spectral atlases 
with complete and consistent optical-NIR coverage are scarce 
for these stellar classes. As argued for the HDF-S, very late 
M-type, methane, T-, and L- dwarfs as well as normal post- 
main-sequence or variable O- to C-rich stars on the asymp- 
totic giant branch are not expected to be as red in the NIR 
as most redshifted galaxy types and/or are very unlikely to 
be present in the small survey areas at high galactic latitude 
covered by the HDF-S and MS 1054-03 fields. As noted by 
iLabbe et al. (2003a), stars are well separated from galaxies in 
the Js-Ks versus I^u-Ks diagram (see also §|8}; our second 
criterion was that the candidates lie on the locus of stars in 
this colour space. We then verified the size of each candidate 
in the WFPC2 /si4 mosaic (non PSF-matched) and excluded 
those that are spatially resolved. The final list includes 84 
stars down to KI°)^q = 24.8 mag. 

7. CATALOGUE PARAMETERS 

The source catalogue is available electronically on the 
FIRES Web site"*. The catalogue format and entries are simi- 
lar to those of the HDF-S catalogue jLabbe et al. I2003al) . and 
described here. 

ID — Running identification number in the catalogue order 
as reported by SExtractor 

x,y — Pixel positions of the objects in the images, based on 
the /Ts-band detection map. The reference is the 4700 x 4700 
frame of the rectified, registered maps at 0."1 pixel"'. 

RA, DEC — Right ascension and declination coordinates 
for equinox J2000.0 (see Eq.|2}- 

/col./, Ccoi,! — Flux measured in the colour aperture (§ I5.2> 
a nd th e associated uncertainty derived from the noise analysis 
(§ 14. 4> in each of the bandpasses i = U,B,V, Veoe ,hi4,Js,H ,Kfi. 
The units are fxJy. 

./iso,!, cTiso,! — Same as above for measurements in the 
isophotal aperture defined by the surface brightness thresh- 
old of fiiK^.As) = 24.8 mag arcsec"^ in the /Tj-band detection 
map (the MAG_ISO computed by SExtractor). 

/tot(^s), o'tot(^s) — j^s-band flux in the total aperture scaled 
by the aperture correction (§ 15. 2> an d the associated uncer- 
tainty from the noise analysis (§ 14. 4> . The units are /iJy. The 
aperture correction is based on the /Q-band PSF profile out- 
side of the total aperture. 

/auto(-Ks), o'auto(-S's) — Same as above for measurements in 
the Kron-like autoscaling aperture based on the light distri- 
bution in the /Tj-band detection map (the MAG_AUTO com- 
puted by SExtractor). 

fap.i.d, fap.i.d — Flux in the bandpass / and t he a ssoci- 
ated uncertainty derived from the noise analysis (§ I4.4> mea- 
sured in each of the 30 circular apertures with diameter d = 
7,8,9, 10,..., 29, 30,35,40,45, 50,55, 60 pixel of the rectified 
maps (0."1 pixel"'). The units are iiiy. 

apco l — Type of aperture used for measuring the colour flux 
(§ \5.2l . encoded as follows. 1: minimum circular aperture 
considered with diameter d=\ ."0 used whenever the isopho- 
tal aperture is smaller 2: maximum circular aperture con- 
sidered with diameter d = 2."0 used whenever the isophotal 
aperture is larger. 3: isophotal aperture defined by the sur- 
face brightness threshold of /i(A's.AB) = 24.8 mag arcsec"^ in 

' * http ://w w w. strw. leidenuniv. nl/~fires 
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the /Ts-band detection map. 4: circular aperture with diameter 
d = -y^Aiso/Tr/i where i = 1.4 is the aperture reduction factor 
applied for blended sources. 

aptot — Type of aperture used for measuring the total K^- 
band flux ('§ I5.2> . encoded as follows. 1: Kron-like autoscal- 
ing aperture based on the light distribution in the K^-bwd de- 
tection map. 2: Colour aperture used in the case of blended 
sources. 

apcorr — Scaling factor corresponding to the aperture cor- 
rection applied to the K^-hwd flux measured in the total aper- 
ture in order to estimate the total integrated flux. 

'"col, '"tot — Circularized radii r = \/A/tt from the area of the 
apertures used to measure the colour and total fluxes. Units 
are pixels in the rectified images (0."1 pixel"'). 

Aiso, Aauto — Area of the aperture enclosed within the 
M(^s.ab) = 24.8 mag arcsec"^ isophote, and of the elliptical 
Kron-like autoscaling aperture (irab, where a and b are the 
semi-major and semi-minor axis lengths). Units are squared 
pixels in the rectified images (0"1 pixel"'). 

FWHM(/:j, FWHM(/8i4) — FWHM of sources from the 
/iTs-band detection map assuming a Gaussian profile, and from 
the non PSF-matched WFPC2 /g^ mosaic. The latter was 
obtained by running SExtractor on the I^u map, and cross- 
correlating the /814- and /STj-band selected catalogues. Units 
are pixels in the rectified images (0"1 pixel"'). 

Wi — Average of the effective weight in the bandpass / 
within the point-source aperture for the PSF-matched set of 
maps (diameter of c/ = 1.5 x FWHM = 1''0). For the optical 
FORS 1 and WFPC2 data, the effective weight is simply pro- 
portional to the total exposure time at each pixel while for the 
NIR ISAAC mosaics, it corresponds to the inverse variance of 
the effective background noise resulting from the integration 
time and the s ystem ati cs fr om the raw data and the reduction 
procedure (§§ l3.1l and l3.2L The weights are normalized to a 
maximum of unity over the maps. 

Wf{Js,H,Ks) — Same as above for the weights in each of 
die fields F = 1,2,3,4 of the NIR ISAAC data separately. 
These weights are simply proportional to the total exposure 
time per pixel, and normalized to a maximum of unity within 
each field map. 

bias, blended, star — Flags indicating the following if set 
to 1 (and otherwise set to 0). Bias is set if the Kron-like au- 
toscaling aperture APER(AUTO) contains more than 10% of 
bad pixels or if the flux measurement within this aperture is 
affected by neighbouring sources as determined from the light 
distribution in the /Q-band detection map. Blended is set if the 
source overlaps with adjacent objects based on the /Tj-band 
detection map. Star is set if the star identification criteria are 
satisfied (SED better fit by a stellar template than a combina- 
tion of galaxy templates, colours lying on the stellar locus in 
the Js-Ks versus diagram, spatially unresolved light 

profile in the non PSF-matched WFPC2 mosaics; §|6}. 

Zph, crup.iow(Zph) — Photometric redshift and associated 68% 
confidence intervals (zph-ffL; §|6}. 

8. ANALYSIS 

8.1. Completeness in the Ks band 

We derived the completeness limits in the /Ts band from 
simulations in which we added point sources to the mosaic 
and analyzed the recovery fraction as a function of input 
magnitude. We performed the simulations on the rectified 
(0"1 pixel"') and non PSF-matched (FWHM = 0."52) image 
used for the source detection. We treated the four ISAAC 



fields separately because of their different effective depths. 
We considered only the deepest part of each field with nearly 
uniform image quality, where the effective weights exceed 
95% of the maximum within the given field (i.e., wp > 0.95). 
We used the PSF profile (created by avera ging the images 
of bright, isolated, unsaturated stars, see § 14. 2> to generate 
artificial point sources with total /Ts.AB-band magnitudes be- 
tween Kill and 28. For realistic simulations, we distributed 
the sources according to the power-law d^V / dK^ cx . We 
adopted a = 0.20 as determined from the raw number counts 
at 21.4 < < 23.4 mag, where S/N > 10 and where in- 
completeness does not yet play a significant role. 

For each field, we added 30000 simulated point sources to 
the image, drawn at random from the power-law distribution 
and placed at random locations. To avoid overlap among the 
simulated sources, we ran the simulations by adding only 25 
of them at a time to the image. For each realization, we ex- 
tracted the sources as described in § 15. II except for the appli- 
cation of a fainter surface brightness detection threshold of 
M(^s,ab) = 25.3 mag arcsec"- (w 3cr of the background noise 
in the two intermediate-depth fields) to ensure reliable deriva- 
tion of the 50% completeness limits. We considered a source 
as recovered if its measured position lay within 3.5 pixels of 
its input position, although experimentation showed that the 
results are little sensitive to the exact choice of this criterion. 
We accounted for blending with real sources following Eq. |6l 
when determining the total /Tj-band magnitudes of the simu- 
lated sources. We did not apply any particular treatment for 
simulated sources that may have been lost on or falsely re- 
covered as real sources ("unmasked case"). To quantify this 
effect, we repeated the analysis excluding all regions within 
the isophotal areas of real sources ("masked case"). 

Figure |9] shows the completeness curves for each field and 
the output versus input total A's-band magnitudes for the un- 
masked case. Table |5] gives the 90% and 50% complete- 
ness limits for both the unmasked and masked cases. For the 
unmasked case, the 90% completeness levels vary by up to 
0.8 mag between the fields, with a mean of /^^"^g w 24. 1 mag. 
The 50% completeness levels are essentially the same for the 
four fields at Kl% = 25.30±0.01 mag (and close to the 3a de- 
tection limits). The shape and slope of the completeness curve 
for the shallowest field F2 at faint magnitudes is distinctly dif- 
ferent from that of the other three fields, which reach fainter 
and more similar depths. This is caused by a more important 
contribution from spurious detections in F2 because of higher 
amplitude noise fluctuations. 

Not only the depth of the images but also the crowding of 
real sources influences the derived limits, through loss and 
confusion of simulated sources overlapping with real objects. 
This effect is significant in the MS 1054-03 field compared to 
the HDF-S because of the presence of the cluster. We found 
this alters mainly the absolute level and slope of the complete- 
ness curves above 90%, and therefore the inferred 90% lim- 
its. For comparison, in the masked case the bright end of the 
curves is flatter, and the 90% completeness levels increase to 
Ri 24.7 mag and reflect more closely the relative depth be- 
tween the fi elds as derived from the background noise analy- 
sis in § 14.41 In contrast, the steepest part of the curves is more 
robust against loss and confusion of simulated sources with 
real sources, with derived 50% limits increasing marginally 
to 25.34 ±0.02 mag. 

Our assumption of point-like simulated sources optimizes 
the detection efficiency, which is further enhanced by our fil- 
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tering of the detection map with a PSF-like kernel (see § 15. U . 
Consequently, the completeness levels derived above consti- 
tute upper limits. In reality, extended sources would have 
brighter completeness levels that would also depend on the 
actual sizes and light profiles. An analysis as a func tion of 
source size and morphology is presented by Truiill oet al. I 
(^005). In correcting the raw counts below, we followed the 
simple conservative approach based on the point-source simu- 
lations presented above, and used the results for the unmasked 
case. 

8.2. Ks-band number counts 

We computed the raw /iTs-band source counts over the same 
area of the mosaic as considered for the simulations presented 
above (where wp > 0.95 for each individual field). We ex- 
cluded stars, used 0.5 mag wide bins, and extended the counts 
to the 50% completeness level (/r'°^g w 25.3 mag). In cor- 
recting the raw counts for incompleteness and spurious de- 
tections, we accounted for the non-uniform depth and noise 
properties over the mosaic as follows. We first treated each 
field separately, computing the raw counts, applying the cor- 
responding corrections, and weighting by the effective area 
where the depth equals the faint edge of each magnitude bin. 
We then combined the corrected counts from the four fields 
and normalized by the total effective area over the mosaic as 
a function of depth to obtain the surface densities. We derived 
the effective area to a g iven depth based on the mosaic's effec- 
tive weight map (Wgff; § 13. U . which accounts for the variations 
in absolute noise levels over the entire field of view. 

This procedure yields the most reliable results for our 
MS 1054-03 data, especially because of the significant off- 
set in central depth of F2 compared to the other fields. Ac- 
counting for the effective area as a function of depth has 
little impact within each field because we only considered 
the deepest regions within each of them. However, it be- 
comes important when combining the fields' counts. In par- 
ticular, at the faintest magnitude bin considered (centered at 
^s°AB - 25.0 mag), the raw counts as well as the point-source 
simulations in F2 suffer from a larger contamination by spuri- 
ous detections so that the derived corrections are less reliable. 
Our treating the fields separately and accounting for the effec- 
tive area per magnitude interval avoids a strong bias from F2 
at the faint end. 

Figure[roishows the differential /Tj -band galaxy counts for 
the MS 1054-03 field. The raw counts are plotted along with 
the counts corrected for incompleteness only and for both in- 
completeness and spurious detections. The incompleteness 
correction is simply the inverse of the detection efficiency of 
simulated sources given by the completeness curves. The cor- 
rection for incompleteness and spurious detections is taken as 
the ratio of the number of simulated sources per input mag- 
nitude bin to the number of detections per output magnitude 
bin. Since spurious sources contribute to the number of de- 
tections (as for the raw counts), the full correction is effec- 
tively smaller than that for incompleteness only. As the figure 
shows, little correction is needed out to the faintest magni- 
tude, where the incompleteness and full corrections amount 
to 30% and 10%, respectively. 

The presence of the z = 0.83 cluster has minimal effect on 
the counts at the faint end. The average magnification fac- 
tor from gravitational lensing over the entire mosaic is small 
and ranges from 1 .05 to 1 .25 for background source redshifts 
between z = 1 and 4 (H. Hoekstra, priv. comm.; based on 
lHoekstraetan i2000). Moreover, the flux and area magni- 



fication cancel out to first order by shifting the counts in the 
surface densities-magnitude plane along a path that is roughly 
parallel to the observed distribution (as indicated in Figure 
for z = 2.5 with average magnification factor of 1.20). At 
brighter magnitudes, the cluster galaxies produce the bump 
around KI°)^q w 20.5 mag (the approximate cluster contribu- 
tion is illustrated in Figure [TO] taken as the counts of all 
sources with 0.6 < z < 1.0, which are dominated by the cluster 
galaxies; see Figure |8}. Accurate corrections for the lensing 
effects and direct contribution to the counts of the cluster are 
beyond the scope of this paper (and would depend, e.g., on 
the luminosity function of the background sources, which is 
unknown, and on a reliable and complete census of the cluster 
members). 

Figure ^3 compares the FIRES counts with those of se- 
lected deep surveys from the literature. In the range « 
22-25 mag (Kl%^„^ w 20-23 mag), the MS 1054-03 field 
counts follow a dlog(A^)/dmag relation with logarithmic 
slope a = 0.20 (it varies from 0.19 to 0.21 between the raw 
and incompleteness-corrected cases). This is comparable to 
the faint-end slope derived from our w 0.7 mag deeper Ks- 
band map of HDF-S ( Labbe et al. 2003a). Overall, the counts 
in the two FIRES fields agree quite well out to the 50% com- 
pleteness limit of the shallower MS 1054-03 field, with the 
exception of the bump due to the MS 1054-03 cluster itself. 
At the faint end, the good agreement suggests that faint cluster 
galaxies probably do not significantly affect the count slope 
in the MS 1054-03 field. The w 5 times wider area of the 
MS 1054- 03 field also leads to more robust statistics com- 
pared to the deeper but smaller HDF-S. Compared to other 
TiT-band surveys, the MS 1054-03 field is not the deepest nor 
the widest, but is unique in its combination of depth and area 
surveyed. 

There are appreciably large differences among the various 
counts shown in Figure In particular, the faint-end slope 
of the FIRES fields generally lies at the flatter end of the 
range from the other surveys (w 0.23-0.36 at 22 < K^ ab < 
25-26 mag, or 20 < K^.yeea 5, 23-24 mag; Diorgovski et alj 
1995; Moustakas et al. 1997; Bershadv. Lowenthal. & Kc^ 
1998; Saraccoetal. 1999; Maihai-a et al. 2001). It seems 
possible that the .Bershadv et al. .(,1998) counts may have been 
overestimated. The significance of the differences between 
studies is however difficult to ascertain in the faint regime, 
where model-dependent corrections become important. A 
complex combination of effects is probably at play, includ- 
ing intrinsic field-to-field variations (which can be large since 
all the fields are relatively small, covering a few arcminutes 
squared to a few tens of arcminutes squared) as well as dif- 
ferences in the photometric measurements (e.g., estimates of 
total magnitudes), the filter bandpasses (K, K^, K'), the deriva- 
tions and application of the corrections to the counts (spurious 
sources contaminations, non-uniform depth, size and mor- 
phology dependencies). This makes any comparison between 
surveys far from straightforward. 

8.3. Colour and magnitude distributions 

FiguresEland^lshow the distribution of Aj-band selected 
sources from the MS 1054-03 field in various Is\4JsHKs 
colour-magnitude and colour-colour diagrams. The sources 
are selected to have a S/N > 10 on their /Q-band colour flux, 
and to lie within the well-exposed areas with wp > 0.3 in 
each NIR band and Weft > 0.3 in each optical map (implying 
^s°AB < 24.4 or /:s°vega < 22.5 mag for the resulting sample). 
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All colours re fer to measurements in the colour apertures, as 
described in § 15.21 The colour distributions are unaffected by 
the (achromatic) gravitational lensing of background sources 
from the z = 0.83 cluster The cluster sequence is seen in all 
plots, in the tight relationships at brighter magnitudes. The 
stars identified in § |6l follow a well-defined blue locus in 
Js-Ks and H-K^ versus /Tj"', and in J^-K^ versus /814-As 
and H-K^ colours, well separated from the distributions for 
galaxies. 

As for the HDF-S (Lab be et al. 12 003a). we find large num- 
bers of sources with very red Js.-K^ colours in the MS 1054- 
03 field. In particular, the criterion /s,vega-^s,vega > 2.3 mag 
effectively selects evolved or dust-obscured galaxies at z > 
2, the "Distant Red Galaxie s" (DRGs, Franx et al. 2003; 
Ivan Dokkum et a l. "2003. 2004). Their average rest-frame op- 
tical colours fall within the range covered by normal galax- 
ies locally, unlike the Lyman-break galaxies (LBGs) selected 
in the o ptical as U -dropouts, which are typically much bluer 
dForster Schreiber et al. 2004; see also Papovich et al. 2001 ). 
This indicates that DRGs have on average higher mass-to- 
light ratios than LBGs at similar redshifts. The surface density 
of DRGs is 1.6 ±0.3 and 2.9 ±0.8 arcmin'^ in the MS 1054- 
03 and HDF-S fields, respectively, at Kfl^^^^ < 22.5 mag, and 
their contribution to the integrated rest-frame V-band lumi- 
nosity is ^ 30%. This suggests that DRGs may be a signif- 
icant constituent of the z ^ 2 - 3 universe in terms of stellar 
mass (Franx et al. 2003; Forster Schreiber et al. 2004; Rud- 
nick et al. , in preparation). Field-to-field variations are likely 
large for DRGs, as suggested by th e high clustering ampli- 
tude derived for the HDF-S sample (Daddi et al. 2003) and 
by the fact that DRGs lie at the high-mass end of galaxies at 
z > 2 (Forster Schreiber et al. 2004). These variations may 
explain why such red J-K selected galaxies have not been 
recognized as an important high redshift population in previ- 
ous NIR surveys; those deep enough to detect very red and 
faint high redshift sources were limited to very small areas. 

Another class of red objects, the Extremely Red Objects 
(EROs) at 1 < z < 2 selected by their very red I - K or: 
R-K colours, ha s been well-studied in the past decade (see 
iMcCarthvl 12004. for a review). Figure [T51 shows that at 
^s°vega < 22.5 mag in the MS 1054-03 field, there is some 
overlap between samples selected with the DRG criterion 
-/s.Vega - ^s,vega > 2.3 and with the FRO criterion /gu^vega- 

■^^s.Vega > 4. About 60% of 7s,Vega--^^s,Vega > 2.3 objects 

also have /gu.vega-^s.vega > 4 colours. However, w 30% only 
of /8i4.vega-^s,vega > 4 objccts meet the DRG criterion. This 
reflects primarily a redshift effect: the DRG criterion is effi- 
cient at isolating evolved or very dusty systems at z > 2, as 
it is intended to, while the FRO criterion is more sensitive to 
galaxies at lower redshifts, down to z ^ 1 . Figure[T4lcompares 
the redshift distributions of the samples in the MS 1054-03 
field. 

9. SUMMARY 

We have presented the observations, reduction, source de- 
tection and photometry of the field around the MS 1054-03 
cluster observed as part of FIRES. We have described the 
main properties of the final images and of the /Tj-selected 
source catalogue. The MS 1054-03 field complements our 
HDF-S data, covering an area w 5 times wider with very sim- 
ilar image quality in the NIR and reaching limiting total mag- 



nitudes for point sources about 0.7 mag brighter. We have 
compared the photometric redshifts with the extensive set of 
spectroscopic redshifts available for the MS 1054-03 field 
and found very good agreement with (|zsp-Zph|/(l + Zsp)) = 
0.078. The MS 1054-03 field confirms our initial findings 
from the HDF-S, with large numbers of red NlR-selected 
high redshift objects, and substantially increases the statisti- 
cal significance of the samples, in particular for the Distant 
Red Galaxies at z > 2 selected by their /s,vega-^s,vega > 2.3 
colours. 

This paper completes the initial phase of the FIRES survey. 
Figure^lplaces FIRES in a broader context by comparing it 
in the area-depth plane to other published deep NIR surveys. 
The HDF-S data constitute the deepest ground-based imaging 
at NIR wavelengths to date. In J and H band, it is only sur- 
passed in depth by space-based imaging with HST NICMOS 
while in ^T-band it represents the currently deepest observa- 
tions even from space. The MS 1054-03 field occupies a 
unique position in terms of area-depth combination, being the 
widest for its depth or conversely the deepest for its size. 

We have exploited our deep optical-NlR imaging of the 
HDF-S and MS 1054-03 fields in several papers to investi- 
gate the properties of the ^Tj-selected galaxies and various sub- 
samples, focusing on their SEDs, their rest-frame optical lu- 
minosities, colours, sizes, their stellar masses, their clustering, 
and the derived redshift evolution of both their integrated light 
and stellar mass (Franx etal. 2003; Rudnicketal. 200^ 
iDaddi et al. 2003; Labbe et a l. i i2003bL .van Dok kum et aU 
2003, 2004; Forster Schreib er et al. I 12004 iTruiillo et all 
I2003i (2005). The results of this initial phase have guided 
the directions of our current and planned follow-up programs. 
The scientific inferences from the HDF-S and MS 1054-03 
fields remain yet limited by the effects of field-to-field vari- 
ations, and the size and low spectroscopic completeness of 
the high-redshift samples in particular We have thus ac- 
tively undertaken the addition of more disjoint fields and ex- 
tensive spectroscopy. Current and planned follow-ups also 
include near-infrared integral field spectroscopy, extension 
of the imaging to other wavebands (especially in the mid- 
IR with Spitzer), and HST ACS and NICMOS high angu- 
lar resolution imaging to constrain the nature and detailed 
spectral, dynamical, and morphological properties of the 
galaxy sam ples (e.g ., Rubin et al. ..2004^ . Labbe et al. , .200^ 
van Dokku m etaLil20 05; Toft et al. Il2005t K. K. Knudsen et 
al. , in preparation). With its unique data set, FIRES has been 
influential to our understanding of galaxy evolution, and will 
continue to be so as it is extended by its various follow-up 
programs. It sets an important basis in the framework of ma- 
jor current and planned surveys such as, for instance, GOODS 
and COSMOS. 
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Colour cutoffs to select EROs vary somewhat in the literature depending 
on the filter bandpasses involved and the scientific purposes of the authors. 
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Fig. 1 . — Fields of view of MS 1054-03 covered by the different instruments. The full I SAAC field of view is outlined by the thick solid contour around 
the A's-band mosaic. The mosaic is normalized by the squared root of the weight map (§ 13. H and displayed with a Unear intensity scaling. The thin solid lines 
delineate the four individual ISAAC fields labeled Fl to F4, and show the regions of overlap. The dashed and dotted lines indicate the outer limits of the WFPC2 
mosaic and the FORSl field of view, respectively. The cluster is apparent as the east- west elongated S-shaped structure straddling fields F2 and F3. 



FIRES: the MS 1054 - 03 field data 17 



T ' ' ■ I ' ' ■ I ' ' ■ I ' ' ■ r 




FWHM [arcsec] 



Fig. 2. — Distributions of the seeing FWHM for the NIR ISAAC observations. The seeing was measured by fitting Moffat profiles to bright, isolated, 
unsaturated stars in the MS 1054 - 03 individual reduced frames. Histograms are shown for each band separately: Js {empty histogram with thick grey line), H 
{hatched histogram), and Ks {grey-filled histogram). 
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Fig. 3. — Selected observational parameters for the NIR ISAAC data set. 7s-, H-, and A's-band data are plotted with Mack triangles, grey crosses, and black 
cirlces, respectively, (a) Relative instrumental source count rates derived from bright stars in the MS 1054-03 Held versus modified Julian date (in days). The 
arrows at the bottom of the plot indicate the February 2000 and March 2001 VLT Antu primary mirror realuminization, separating the periods PI, P2, and P3 of 
data acquisition. The count rates were integrated in a 6"-diameter aperture and normalized to the median values over period P3. (h) Relative stars count rates 
versus airmass. (c) Background levels in the individual raw frames versus modified Julian date. The background corresponds to the median count rates in the 
raw frames. Periods PI to P3 are identified at the bottom of the plot, (d) Background levels versus airmass. (e) Background levels versus Universal Time (in 
hours). (/■) Seeing FWHM of individual frames versus airmass. In all plots some OBs depart from the general trends, reflecting particular conditions during the 
observations. 
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Fig. 4. — Effects on the PSF of the combination and rectification of the ISAAC OB images. The FWHM of the PSF measured in the rectified combined OB 
images is compared to the median of the FWHMs in the individual non-rectified frames within each OB. The angular resoluti on is preserved to within ±5% at 
FWHM > Of'5 and degrades progressively towards lower FWHM by up to 25% at ft: Of'3, at the Hmit of undersampling (see ii l3.H . The FWHM was measured 
by fitting Moffat profiles to bright, isolated, unsaturated stars. Data are shown separately for the three NIR bands: J^: black triangles; H: grey crosses; Ks; black 
circles. 
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Fig. 5. — Three-colour composite image of the central part of the MS 1054-03 field. The image is constructed from the WFPC2 /814 map (coded in blue), 
and the ISAAC Js and Ks maps (coded in green and red, respectively). All three maps are convolved to a common angular resolution of FWHM = 0"69 (the 
PSF of the l/-band map, which has the poorest resolution among the full data set). The field of view shown is 4' X 4'; North is up and East is to the left. There 
is a rich variety of colours among the sources, especially for the fainter ones. Several of the red objects have spectroscopic or photometric redshift z > 2, and 
are evolved, dust-obscured, massive galaxies (Forster Schreiber et al. 2004). Most are extremely faint in the observer's optical bands and would be missed in 
optically-selected high-z samples. The WFPC2 mosaic does not cover the top left and right corners, which causes the colours of many sources to appear yellow 
in these regions. 
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Fig. 6. — Empirical scaling relation of the baclfground noise rms of the MS 1054-03 field data set as a function of linear size N = \/A for apertures of 
area A. The plots show the relations for the four fields of the ISAAC A's-band mosaic used for the photometric measurements (rectified and resampled to the 
WFPC2 scale of 0" 1 pixel"' and smoothed to match the PSF of the FORSl f/-band data with FWHM = Of'69). The results are qualitatively similar for all other 
bands, (a) Measured rms variations of the background noise for different aperture sizes in each of the four ISAAC fields (Fl: filled dots; F2: open dots; F3 : 
open squares; F4 : filled triangles). The solid lines indicate the best-fits to tlie measured relations following Eq.|3] The dashed lines show the linear predictions 
from the measured pixel-to-pixel rms if the noise were purely uncorrelated and Gaussian, (b), (c), (d), and (e) Distributions of count rates measured in random 
apertures within empty areas of the four ISAAC fields maps (histograms), which are used to determine the background noise rms for each aperture size. The 
histograms are plotted for three selected aperture sizes with N = 8.86, 13.3, and 17.7 pixel, coiTesponding to circular apertures of diameters = 1", lf'5, and 2". 
The background noise rms is derived from the best-fit Gaussian profile to the distributions (black dashed lines). 
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Fig. 7. — Comparison of various metliods to estimate total Ki-band magnitudes. Tlie plots show the difference between the various estimates and the expected 
total magnitude for point sources, given by the magnitude in a If'O-diameter circular aperture Ks(l") scaled by the total aperture coirection of 0.58 mag (see 
§15^. The values of the abcsissa are our adopted total ATs-band magnitudes. Stars are distinguished from galaxies {grey stars and black dots, respectively, see 
§161. The error bars correspond to Icr uncertainties. For point sources, the best A^J°' estimates should be equal to ^Gi(l") — 0.58 over the entire range in brightness 
(horizontal black dashed lines), (a) Isophotal magnitudes K^"" strongly underestimate the total magnitudes Kl"' at the faint end where most sources are point-like 
and where extended emission lies below the isophotal surface brightness threshold for detection and photometry, (b) Kron-like autoscaling magnitudes K'"^" still 
underestimate the of the fainter objects by up to 0.58 mag, reflecting the emission missed from the PSF wings in the smaller apertures, (c) The turnover at 
faint magnitudes is nearly absent when using the Kl°' as measured in this paper 
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Fig. 8. — Photometric and spectroscopic redshifts in the MS 1054-03 field, (a) Direct comparison between the photometric redshifts Zph and spectroscopic 
redshifts Zsp for all Xj-band detected sources at ArJ°^g < 24.8 mag for which a Zsp measurement is available. The vertical error bars indicate the 68% confidence 
intervals on Zph derived from Monte-Carlo simulations (see §|^. In the ideal case, Zph = Zsp {gi~ey solid line), (b) Residuals dz/(l +Zsp) = (Zph — Zsp)/(1 +Zsp) 
as a function of spectroscopic redshift. The accuracy of the photometric redshifts quantified by the mean {|dz|/(l +Zsp)} is 0.078 (dotted lines) about the ideal 
value of (solid line), (c) Distributions of photometric redshifts for all A's-band detected sources at K'°j^^ < 24.8 mag (hollow histogram) and of spectroscopic 
redshifts for the subset of objects for which they are available (grey filled histogram). In all plots, stars have been excluded. 



24 Forster Schreiber et al. 




Fig. 9. — Results of the completeness simulations for point sources in the Ks band for the MS 1054-03 field. For the plots shown here, the analysis was 
carried out considering the entire area in each of the ISAAC fields with effective weight wp > 0.95 including also the isophotal areas of real sources, (a) 
Completeness curves giving the fraction of recovered simulated sources as a function of input total Ks-band magnitude. The results for each ISAAC field are 

plotted separately (black line and filled circles: Fl; grey line and open circles: F2; grey line and open squares: F3; black line and filled triangles: F4). The 90% 
and 50% completeness limits are indicated {upper and lower dashed lines), (b) Observed versus input total -Ks-band magnitude for all artificial point sources in 
the simulations in all four ISAAC fields. The ideal relation where the observed magnitudes equal the input magnitudes is indicated (grey dashed line). 
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Fig. 10. — Differential Xs-band gala xy co unts in the MS 1054 — 03 field. The counts are given as a function of the total Xs-band AB magnitude (with total 
magnitudes computed as described in § 15.21 . The raw counts {grey open circles) are summed up in bins 0.5 mag wide (plotted at the central bin values) and 
over the total area where each ISAAC field has an effective weight wp exceeding 95% of its maximum. The counts were corrected for incompleteness (grey 
crosses) and for both incompleteness and spurious detections (black filled circles) using point-source simulations. Because the four fields have different noise 
properties, the corrections were derived and applied for each field separately, and further account for the variations of effective area as a function of depth. The 
corrections are small except in the faintest bin near the 50% completeness or 3cr detection limit (dashed vertical line). Correction for incompleteness only 
overestimates the counts at the faint end where the contribution from positive noise peaks becomes important. The counts at ^°^g ~ 22 — 25 mag (Ki 20-23.5 
in Vega magnitudes) have a logarithmic slope of 0.20 (grey solid line). The error bars represent Poisson uncertainties. Galaxies from the z = 0.83 cluster 
produce the bump seen around K'°j^^ Ri 20.5 mag, as indicated by the number counts of sources around the cluster redshift (dotted line). Gravitational lensing by 
the cluster affects little the surface densities: the effects are small on average over the area surveyed and act along a di rection roughly parallel to the counts curve 
as illustrated for the mean magnification by 20% for z = 2.5 background source plane (arrow at bottom left; see 8 18.21 . 



M 
D 
Xi 



105 



QO 
CO 



I ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 

• MS1054-03 field 
■ HDF-S 

* Maihara et al. 8001 
J- Saracco et al. 1999 
X Bershady et al. 199B 
+ Moustakas et al. 1997 
A Djorgovski et al. 1995 



10^ 



00 

o 



103 



16 17 18 19 20 21 22 23 24 
ifvega [mag] 



Fig. 1 1. — FIRES A's-band galaxy counts compared to published counts in other deep A^-band fields. All data are corrected counts, and are here plotted as 
a function of magnitudes in the Vega photometric system. The legend gives the symbols used for different references. All counts are plotted up to their 3(7 
detection or 50% completeness limits, except for the Saracco et al. ( 1999) data which extend to a S/N 5. Error bars are plotted for the FIRES counts and 
represent Poisson uncertainties. The 5 times larger area of the MS 1054 — 03 field leads to more robust statistics than for the deeper but smaller HDF-S up 
to /Tvcga ~ 23 mag, a depth only about 1 mag brighter than current deep AT-band surveys. The cluster bump around A^veea = 18-5 mag is conspicuous for the 
MS 1 054 - 03 field. The counts slope at AVcsa > 20 mag in both FIRES fields Ues generally at the flatter end of results reported in the literature. The differences 
between various surveys are probably due in part to field-to-field variations, and different methods used for the photometric measurements and for the derivation 
of the corrections to the raw counts. 
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Fig. 12. — Colour-magnitude diagrams of S's-band selected objects in the MS 1054-03 field. Magnitudes and colours are on the Vega system. Only sources 
with S/N > 10 on their A^s-band colour flux and a minimum of 30% of the total exposure time in at least one ISAAC field for each NIR band and in the /gi4-band 
map are included. Galaxies at z < 2 and z > 2 as well as stars are shown with different symbols (filled circles, open circles, and open stars, respectively). The 
error bars represent the la photometric uncertainties, and sources with S/N < 2 in their Js-, H-, or /g^-band fluxes are plotted at their 2cr confidence level. The 
dotted line in each panel indicates the limits corresponding to the 3cr background noise level for the point-source aperture in the PSF-matched images used for 
the photometry, (a) Js-Ks colour versus total ^s-band magnitude, (b) H-Ks colour versus total Xj-band magnitude, (c) hu-Ks colour versus total ^s-band 
magnitude. 
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Fig. 13. — Colour-colour diagrams of A's-band selected objects in the MS 1054—03 field. The colours are on the Vega system. Only sources with S/N > 10 on 
their Xj-band colour flux and a minimum of 30% of the total exposure time in at least one ISAAC field for each NIR band and in the /814-band map are included. 
Galaxies at z < 2 and z > 2 as well as stars are shown with different symbols {filled circles, open circles, and open stars, respectively). The error bars represent 
the 1(T photometric uncertainties, and sources with S/N < 2 in their 7s-, H-, or /gM-band flux are plotted at their 2a- confidence level, (a) Js — Ks versus /814 —Ks 
colours. The criterion 7s —Ks > 2.3 mag defined to identify Distant Red Galaxies at z > 2 (DRGs; Franx et al. 2003) and the cutoff Ig[4 — Ks > 4.0 mag to select 
Extremely Red Objects (EROs) are indicated (horizontal and vertical dotted lines, respectively), (b) Js -Ks versus H-Ks colours. The criterion 7s -^^s > 2.3 mag 
to select DRGs at z > 2 is again indicated (horizontal dotted line). 
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Fig. 14. — Redshift distributions of red A's-band selected galaxies in the MS 1054 — 03 field. Only galaxies with S/N > 10 on their S's-band colour flux 
(implying K'°y^^.^ < 22.5 mag) and a minimum of 30% of the total exposure time in at least one ISAAC field for each NIR band and in the /814-band map are 

included. The spectroscopic redshifts were used wherever available. The colour cut /s.Vega ~ ^s, Vega > 2.3 mag is devised to identify evolved galaxies at z > 2 
(Distant Red Galaxies or IDRGs; Franx et al. 2003). The galaxies that satisfy this criterion (filled histogram) all lie at z > 1.8 except for one interloper at z = 1.19, 
with median redshift (z)nicd = 2.36 (2.42 when excluding the interloper). The galaxies that have Iguyega ~^s.Vega > 4 colours of Extremely Red Objects (EROs; 
hatched histogram) typically lie at lower redshifts, with (z)nied = 1.65. There is overlap between the two colour-selected samples (cross-hatched histogram) but 
the /s. Vega ~^s. Vega > 2.3 Criterion selects more efficiently z > 2 galaxies. 
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Fig. 15. — Comparison of published deep NIR imaging surveys in the literature as a function of depth and area, (ci) /-band surveys, including HST NICMOS 
imaging in the Fl lOW bandpass. (6) //-band sureeys, including HST NICMOS imaging in the F160W bandpass, (c) /T-band surveys, including imaging with K, 
Ks, and K' filters. Ground-based and space-based surveys are distinguished by different symbols (filled circles and open circles, respectively). In addition to the 
FIRES HDF-S and MS 1054-03 fields (highhghted by the shaded area), the surveys considered include the following: the NICMOS fields of the HDF-N (HDFN, 
from Thompson et al. 1999 for the deeper but smaller area, and from Dickinson et al. 2003 for the wider and shallower area); the NICMOS field of the Hubble 
Ultra Deep Field (UDF Thompson et al. 2005); the HDF-S parallel NICMOS field (HDFS-NIC-P, Williams et al. 2000); the two UDF parallel NICMOS fields 
(UDF-NIC-P, Bouwens et al. 2005); the Subaru Deep Field (SDF, Maihara et al. 2001); several blank fields imaged with the Keck telescopes (three fie lds in K 
from Dioraovski et al. 1995 and two in / and K fromBershadv et al. 1998); the Chandra Deep Field South as part of GOODS (GOODS-CDFS, Giavalisco et al. I 
|2g04); the ESO New Technology Telescope Deep Field (NTTDF, Saracco et al. 1999); the FORS Deep Field (FDF Heidt et al. 2003); the NICMOS parallel 
imaging database for Cycle 7 of high Galactic latitude fields (NICMOS-Ps, Corbin et al. 2000); the two fields as basis for the K20 survey i Cimatti et al. 2002}); 
the ESO Imaging Survey Deep fields (EIS-Deep, Arnouts et al. 2001); the Calar Alto Deep Imaging Survey (CADIS, Meisenheimei 1998); the Las Campanas 
Infrai-ed Survey (LCIRS, Chen et al. 2002); the Munich Near-Infrared Cluster Survey (MUNICS, Drory et al. 2001). For GOODS-CDFS, we used the average 
depth reported by ESO for the version 1.0 pubhc release of the NIR ISAAC data covering 131 arcmin^ in / and /Ts (shown in black); the expected ai'ea and 
additional //-band data once the survey is completed are also indicated (grey filled circle). As much as possible, the limiting depths refer to total magnitudes, 
or magnitudes in fixed apertures significantly larger than the PSF, for point sources. The FIRES data of HDF-S currently represents the deepest ground-based 
survey in the / and H bands, and the deepest survey in the K band from either ground or space. The wider FIRES MS 1054 — 03 field is the deepest for its size, 
filling the regime between the ultra-deep surveys over smaller areas and significantly larger but yet shallower fields such as GOODS-CDFS. 
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TABLE 1 

Summary of the MS 1054-03 field observations 



Camera 


Filter 


Pointings " 






PSF'' 










(s) 


(arcsec) 


All bands 


FORSl 


U 


1 


30 


16500 


0.69 


FORSl 


B 


1 


18 


7200 


0.57 


FORSl 


V 


1 


6 


1800 


0.65 


WFPC2 


F606W 


6 


12 


6500 


0.21 


WFPC2 


F814W 


6 


12 


6500 


0.22 


ISAAC 


/s 


4 


776 


93120 


0.48 


ISAAC 


H 


4 


739 


87780 


0.46 


ISAAC 




4 


1588 


95280 


0.52 


Near-infrared bands 


ISAAC 


h 


Fl 


205 


24600 


0.44 






F2 


211 


25320 


0.46 






F3 


180 


21600 


0.45 






F4 


180 


21600 


0.46 


ISAAC 


H 


Fl 


196 


23520 


0.46 






F2 


180 


20700 


0.42 






F3 


184 


22080 


0.43 






F4 


179 


21480 


0.42 


ISAAC 




Fl 


473 


28380 


0.46 






F2 


398 


23880 


0.50 






F3 


358 


21480 


0.41 






F4 


359 


21540 


0.39 



" Number of pointing positions. For the ISAAC observations, fields Fl, F2, F3, and F4 correspond 
to the north, east, south, and west quadrants of the combined mosaic, respectively. The properties of 
the entire mosaic and of each field separately are listed. 

^ Number of individual exposures, or frames. Frames rejected because of too poor quality are not 
included (3% and 33% of the total number of exposures taken with ISAAC and FORSl, respectively). 

Total integration time. Data rejected because of too poor quality are not included (3% and 31% of 
the total integration time spent with ISAAC and FORSl, respectively). 

FWHM of the best-fitting Moffat profile to the PSF, determined by averaging the profiles of a 
sample of bright, isolated, unsaturated stars in the combined maps. 



TABLE 2 

Adopted photometric zero points for the MS 1054-03 data set 



Data set Zero point 





Vega mag 


AB mag 


FORSl U 


24.127 


24.828 


FORSl B 


26.746 


26.634 


FORSl V 


27.131 


27.145 


WFPC2 V606 


22.900 


23.020 


WFPC2/8I4 


21.660 


22.090 


ISAAC 7s 


24.970 


25.872 


ISAAC H 


24.846 


26.226 


ISAAC K, 


24.309 


26.169 



Note. — FORSl and ISAAC measurements were performed in a 6"- 
diameter circular aperture. WFPC2 values refer to an "infinite" aperture, 
which is defined as having 1.096 times the flux in an aperture 0"5 in ra- 
dius. The zero points for the FORSl optical bands were derived applying 
mvega = 'Kinstr + ZPvega - AC ■ z + 7c ' Cvcga, where OTvcga is the calibrated magni- 
tude of the standard star, m\nni = -2.5 log(counts rate) is the instrumental mag- 
nitude, z is the airmass, and C is the reference colour for the colour term. The 
extinction coefficients used were n(U) = 0.41, k(B) = 0.21, and k(V) = 0.16, 
in units of magairmass"' . The colour coefficients were ^u-b{U) = 0.10, 
7s-v(6) = -0.09, and 7s-v(V) = 0.05. Extinction an d col our terms were not 
applied for the ISAAC zero points as explained in § 13.11 The WFPC2 zero 
points were directly taken from the WFPC2 manual. 
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TABLE 3 

Background noise in the PSF-matchedMS 1054-03 field images 



Data set ^ 


a" 






1(T background noise limit ind = l"035 




(counts/s) 






(AB mag) 



FORSl U 




1.13 X 10"^ 


1.82 


0.025 


29.01 


FORSl B 




2.25 X 10~^ 


2.31 


0.093 


29.60 


FOP 1 V 




9 24 X 10~^ 


2 30 




xo. / o 


WFPC2 V606 




1.23 X 10^ 


2.01 


0.234 


28.84 


WFPC2 /814 




1.01 X 10^ 


2.19 


0.224 


28.10 


ISAAC 


Fl 


0.00504 


2.57 


0.121 


27.80 




F2 


0.00591 


2.13 


0.140 


27.70 




F3 


0.00531 


1.96 


0.143 


27.87 




F4 


0.00489 


2.65 


0.109 


27.84 


ISAAC H 


Fl 


0.01054 


3.29 


0.101 


27.20 




F2 


0.01192 


2.79 


0.128 


27.14 




F3 


0.00957 


3.06 


0.111 


27.34 




F4 


0.00939 


2.71 


0.136 


27.40 


ISAAC 


Fl 


0.00926 


2.45 


0.128 


27.45 




F2 


0.01323 


2.54 


0.115 


27.07 




F3 


0.00985 


2.76 


0.125 


27.30 




F4 


0.00975 


2.73 


0.143 


27.28 



The images are the rectified smoothed maps at O" 1 pixer' and common FWHM = 0"69 used for 
the photometric measurements. 

^ Pixel-to-pixel rms as measured directly in empty parts of the images, within randomly placed 
circular apertures of 1" -diameter. 

Parameters of the best linear fit to the deviations of the effective backg round rms variations relative 
to the uncorrelated Gaussian scaling as a function of aperture size (§ I4.4> . To estimate the photometric 
uncertainties in a given aperture of linear size A' = sqrt(A) where A is the aperture area, Eq. 0is 
applied with the coefficients listed here and the mean weight within the aperture. 

Effective la background-noise limiting magnitudes of the PSF-matched images in the point- 
source aperture appropriate for the convolved images. The aperture diameter is d = 1.5 x FWHM = 
lf'035 maximizing the S/N ratio in the photometry of point sources. 



TABLE 4 

Background noise in the non PSF-matched MS 1054-03 field images 



Data set " 




Point-source aperture diameter 
(arcsec) 


Icr limiting depth 
(AB mag) 


3(T total limiting magnitude 
(AB mag) 


FORSl U 




1.035 


29.01 


27.24 


FORSl B 




1.035 


29.60 


27.83 


FORSl V 




1.035 


28.78 


27.01 


WFPC2 V606 




0.31 


30.00 


28.28 


WFPC2 /814 




0.33 


29.24 


27.46 


ISAAC 7s 


Fl 


0.72 


28.00 


26.20 




F2 


0.72 


27.77 


25.97 




F3 


0.72 


27.95 


26.15 




F4 


0.72 


27.98 


26.18 


ISAAC H 


Fl 


0.69 


27.39 


25.57 




F2 


0.69 


27.31 


25.49 




F3 


0.69 


27.61 


25.79 




F4 


0.69 


27.63 


25.81 


ISAAC 


Fl 


0.78 


27.61 


25.73 




F2 


0.78 


27.16 


25.28 




F3 


0.78 


27.50 


25.63 




F4 


0.78 


27.54 


25.67 



* The images are the rectified maps at 0" 1 pixer' at their respective intrinsic angular resolution 
(listed in Table0, except for the FORSl data. The B- and V-band images have a similar resolution 
as the y-band image used as reference for the PSF-matching and the differences in the noise analysis 
from the smoothed and unsmoothed FORS 1 sets can be neglected. 

^ Depth of the images as given by the Icr background-noise limiting magnitudes in the point-source 
apertures appropriate for each image, with diameter li = 1 .5 x FWHM. 

Depth of the images as given by the 3a background-noise limiting total magnitudes, i.e., corrected 
for the flux missed outsid e of the aperture based on the growth curve of the PSFs contructed from 
bright stars profiles (see § 15. 2> . 
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TABLE 5 

Ks-BANB COMPLETENESS LIMITS FOR THE MS 1054-03 FIELD 



Field" 


All image area 


Masking 


real sources 




90% limit 


50% limit 


90% limit 


50% limit 




(AB mag) 


(AB mag) 


(AB mag) 


(AB mag) 


Fl 


24.33 


25.30 


24.76 


25.34 


F2 


23.91 


25.31 


24.60 


25.36 


F3 


23.73 


25.29 


24.64 


25.32 


F4 


24.51 


25.31 


24.73 


25.34 



" The A's-band image used for the completeness simulations is 
the rectified, non PSF-matched image at 0"l pixer' and with 
FWHM = 0."52 used for the source detection. 

Completeness limits derived in simulations where the entire 
area in each field with effective weight wf > 0.95 is considered, 
including the regions within the isophotal areas of real sources. 

Completeness limits derived in simulations where the area in 
each field with effective weight wf > 0.95 and outside of the 
isophotal areas of real sources is considered. This reduces the 
effects of loss and blending of artificial sources coinciding spa- 
tially with real sources. The impact is significant for the 90% 
limits, which shift towards fainter magnitudes. 



